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SUMMARY 


OBJECTIVE 

The  objective  is  to  illustrate  the  advantages  arising  from  a  rational  transition 
process  when  developing  an  electrochemical  power  source  designed  to  operate  at  its 
limit.  Here,  we  summarize  our  approach  and  results  based  on  such  a  transition 
philosophy.  Although  we  selected  the  LilS0Cl2  couple  for  this  purpose,  such  an 
approach  is  applicable  to  other  systems  as  well  and  is  an  example  of  current  trends 
in  the  development  of  electrochemical  cells. 

RESULTS  AND  CONCLUSIONS 

In  the  early  1980s,  the  Lij SOCI2  system  was  a  candidate  power  source  for  propul¬ 
sion.  A  major  effort  was  undertaken  to  develop  this  power  source.  To  maximize  the 
power  density  of  this  power  source,  two  design  approaches  were  suggested:  the  con¬ 
struction  of  a  module  containing  in-series  connected  thin  cells,  as  one  choice,  and  the 
employment  of  a  flowing  electrolyte,  as  the  other.  To  improve  the  performance  in  each 
of  the  proposed  designs,  modifications  of  both  electrolyte  composition  and  cathode 
structure  and  composition  were  undertaken,  with  mixed  results.  Optimization  of  this 
system  was  not  achieved  and  the  programs  developing,  this  source  were  terminated. 

The  Lij  SOCI2  system  is  an  example  of  technology  prematurely  transitioned  from 
directed  research  to  battery  development.  In  this  communication  we  outline  a  rational 
approach  to  battery  development.  The  Li/SOCl2  technology  has  come  a  long  way 
since  the  initial  recognition  of  its  value  as  a  battery  chemistry.  Much  of  this  progress 
has  resulted  from  a  focused  effort  to  pursue  a  logical  progression  of  transitions  through 
the  R&D  cycle.  These  efforts  have  elucidated  the  chemistry  and  electrochemistry  of 
the  system,  the  optimal  design  features  of  cells  and  modules,  the  assembly  of  these 
cells  or  modules  into  batteries,  and  the  resultant  characteristics  of  these  assemblies. 
Although  the  Li/SOCl2  system  is  used  as  an  example,  other  battery  technology 
developments  would  benefit  from  such  an  approach. 
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1.  LITHIUM  BATTERY  TECHNOLOGY 


1.1.  INTRODUCTORY  NOTES 


Today,  electrical  energy  is  indispensable.  In  the  broad  concept  of  energy  man¬ 
agement,  batteries  occupy  a  prominent  role  due  to  their  very  high  efficiencies  of 
conversion  and  their  ability  to  act  as  portable  energy  storers.  A  battery  is  probably 
the  most  versatile  electrical  energy  source  both  in  its  performance  and  construction. 
It  covers  a  wide  range  of  power  outputs;  from  10“^  W  needed  to  power  a  wrist  watch 
to  ca  10^  W  required  for  propulsion  applications.  It  can  be  constructed  in  various 
sizes:  from  a  button  size  for  hearing  aids,  to  large  installations  for  load  leveling.  Fur¬ 
thermore,  it  can  deliver  this  power  in  millisecond  pulses  or  continuously  for  several 
years. 


1.2.  EMERGENCE  OF  LITHIUM  BATTERY  TECHNOLOGY 


The  past  few  decades  have  seen  a  rapid  growth  in  battery  technology.  In  addition 
to  time-honored  and  proven  lead  acid  storage  battery,  new  types  of  electrochemical 
power  sources  have  appeared.  Among  them,  perhaps  faster  than  expected,  are  lithium 
{Li)  batteries,  which  have  emerged  as  a  new  class  of  energy-dense  cells.  The  emer¬ 
gence  of  Li  as  an  anode  material  can  be  traced  to  early  research  at  the  University  of 
California,  Berkeley,  where  Harris  (1958)  showed  that  Li  could  be  electrodeposited 
from  solutions  of  lithium  salts  dissolved  in  propylene  carbonate.  The  energy-dense 
LifSOCl2  system,  first  disclosed  in  1973  (Auborn  et  ah,  1973),  was  borne  unex¬ 
pectedly  as  a  “spin-off”  from  fundamental  research  on  new  concepts  in  liquid  lasers, 
a  program  sponsored  by  the  Office  of  the  Chief  of  Naval  Research.  Following  this 
disclosure,  considerable  research,  including  applied  research  and  exploratory  develop¬ 
ment  led  to  an  early  transition  to  engineering  development  about  1980.  Somewhat 
later,  a  full-scale  effort  to  construct  a  multimodule,  high-discharge-rate  battery  was 
underway. 


The  LA/ SOCI2  system  is  an  example  of  rapidly  expanding  electrochemical  power 
source  technology.  Since  its  discovery,  the  battery  has  been  developed  for  use  in  a 
number  of  diver.se  specialty  applications.  Its  main  feature  is  that  it  has  the  highest 
energy  density  of  any,  as  yet,  available.  Recent  advances  in  the  lithium  battery  tech- 
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nology,  including  Lij SOC12-,  have  been  treated  in  a  number  of  monographs  discussing 
various  aspects  of  this  technology  (Salkind,  1987;  Dey,  1987;  Venkatasetty,  1984;  Ga- 
bano,  1983).  While  much  progress  has  been  made,  a  number  of  issues  remain  and 
must  be  addressed  before  the  system  is  ready  for  full-scale  applications.  In  our  opin¬ 
ion,  the  activities  associated  with  the  development  process  itself,  i.e.,  the  transition 
from  research  concepts  to  the  construction  of  engineering  models,  have  not  been  suf¬ 
ficiently  emphasized.  Consequently,  we  focus  on  those  aspects  that  were,  and  are,  of 
interest  to  us,  viz.,  application  of  engineering  principles  when  evaluating  and/or  con¬ 
structing  an  electrochemical  converter.  In  fact,  in  an  earlier  communication,  Szpak 
and  Smith  (1990)  discussed  the  Li/SOCl2  battery  development  program,  identified 
several  problems  and  issues  associated  with  it,  and  proposed  a  systematic  approach 
to  meet  design  objectives. 


In  preparing  this  report,  we  share  our  experience  gained  during  the  development 
process  of  a  specific,  highly  specialized  LifSOCl2  battery,  one  designed  for  high- 
discharge  rates  (Szpak  k.  Driscoll,  1987).  The  need  for  a  systematic  approach  is 
evident  considering  the  energy  density  of  advanced  lithium  systems  is  on  the  order 
of  30-50  percent  of  common  explosives.  This  high-energy  density  simultaneously  of¬ 
fers  great  flexibility  to  the  designer  and  presents  potential  danger  to  the  user.  The 
issues  addressed  here  range  from  research  topics  dealing  with  electrocatalysis,  kinet¬ 
ics  of  electroredu'^iion  of  SOCI2,  physico-chemical  properties  of  electrolytes,  growth 
habits  of  LiCl  crystallizing  within  the  confines  of  the  porous  structure  of  the  positive 
plate  to  engineering  aspects  such  as  thermal  management,  scale-up  problems,  inter¬ 
cell  currents,  dimensional  tolerances,  and  quality  control  measures.  All  these  topics 
are  developed  in  the  context  that  a  battery  is  an  electrochemical  reactor  operating 
in  response  to  imposed  demands  wherein  reactants  are  converted  to  products  in  a 
controlled  manner  with  the  release  of  electrical  energy  and  heat.  When  applying  this 
concept  to  the  Li/SOCl2  battery,  not  only  are  the  overall  chemical  reactions  of  in¬ 
terest,  but  also  the  partial  reactions  taking  place  within  well-defined  sections  of  the 
reaction  space.  An  understanding  of  these  partial  reactions  provide  data  needed  to 
design  a  safe  and  efficient  power  source. 


1.3.  THERMODYNAMIC  CONSIDERATIONS 


A  detailed  discussion  of  the  thermodynamic  structure  of  electrochemical  systems 
was  given  l)y  van  Rysselberglie  (1966)  and  that  of  the  Li/SOCl2  by  Godshall  and 
Driscoll  (1984).  In  its  simplest  ca.se,  an  electrochemical  cell  is  a  constrained  system 
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that  reacts  spontaneously  as  the  constraints  are  removed.  The  thermodynamic  prop¬ 
erties  of  the  reactants  and  products  give  a  measure  of  the  energy  available  from  the 
system.  The  enthalpy  of  the  reaction,  AH,  is  the  total  energy  released  while  Gibbs 
free  energy  of  reaction,  AG,  is  the  maximum  amount  of  energy  that  can  be  converted 
to  electrical  energy.  These  quantities  are  related  by 

AG  =  AH  -TAS.  (1) 


Two  additional  thermodynamic  quantities  are  of  interest:  first,  the  thermody¬ 
namic  potential,  Uq  =  —AG/zF,  and  second,  the  thermoneutral  potential, 
Ut  =  —AH/zF.  The  open  circuit  potential,  t4(0),  is  the  potential  measured  be¬ 
tween  the  positive  and  negative  electrodes  under  the  condition  of  no  current  flow,  i.e., 
f/c(0)  =  {U+  —  This  quantity  is  always  positive  and  depends  only  on  the  sys¬ 

tem  chemistry  and  not  on  the  physical  size  of  a  cell.  The  measured  f/c(0)  cannot 
exceed  Uo  if  the  stoichiometry  of  the  reaction  is  properly  written.  The  thermoneutral 
potential,  Ut,  is  of  practical  importance  because  it  provides  the  basis  for  calculating 
•the  efficiencies  of  discharging  cells,  and  is  the  basis  for  the  thermal  management  of 
battery  systems.  The  total  heat  released,  when  the  cell  operates  at  a  potential  very 
close  to  Uc{0),  is  the  thermodynamically  unavailable  energy  for  conversion,  which  is 
given  by  q^nt  =  —TAS.  This  quantity  is  also  expressed  by  qent  =  zF{Ut  —  Uq).  In 
normal  battery  operation,  the  cell  potential  Uc  is  less  than  Uq  and  the  cell  produces 
additional  heat  equal  to  qd  =  zF{Uq  —  Uc).  Therefore,  the  total  heat  released  is  given 

by 

qceii  =  zF{Ut  -  Uc).  (2) 

Conclusions  based  on  applying  eq.  (1)  are  valid  only  when  the  thermodynamic 
functions  are  accurately  known.  To  illustrate,  early  in  the  development  period  of  the 
Li/ SOCI2  cell,  it  was  reported  that  AS  >  0  for  the  cell  reaction.  Consequently,  it 
was  thought  that  thermal  management  would  be  no  problem  and  that  the  system 
would  cool  during  the  discharge  (Schlaikjer,  Goebel  &  Marincic,  1979).  This  was 
shown  not  to  be  the  case,  however,  and  heat  is  generated  in  the  course  of  battery 
operation  (Godshall  &  Driscoll,  1984). 

1.4.  WHY  ELECTRICAL  ENERGY  CONVERSION? 


Often  the  performance  of  galvanic  cells  is  compared  with  the  performance  of  other 
systems.  For  such  comparisons  to  be  meaningful,  it  must  start  from  the  same  form 
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of  energy  anc.  terminate  with  the  same  usable  form  (Bronsted,  1955).  To  illustrate, 
we  consider  the  generation  of  electrical  energy  by  the  combustion  of  fossil  fuels  and 
contrast  it  with  the  production  by  galvanic  cells  (Bockris  &  Srinivcisan,  1969).  In 
both  cases,  the  initial  form  is  the  chemical  energy,  —AH.  In  thermal  combustion, 
the  stored  energy  is  converted  into  electrical  energy,  Uc(i)e~,  by  the  following  set  of 
consecutive  transformations: 

-  AH  -^Qt^  W{pV)  Ucii)e-.  (3) 


The  intermediate  step,  W{pV),  is  necessary  because  the  direct  conversion  of  heat 
to  another  form  of  energy  is  impossible.  The  inefficient  step  is  the  Qt  — >  \V{pV) 
conversion  because  thermal  energy  can  only  be  utilized  through  its  transfer  to  an 
auxiliary  system  which,  in  turn,  can  deliver  work,  e.g.,  mechanical  compression  of 
stored  gas  from  which  heat  is  rejected  at  lower  temperature.  The  thermodynamic 
efficiency  of  an  ideal  heat  engine  is  e/i  =  1  —  72/7i,  where  T\  is  the  fluid  intake 
temperature  and  T2  its  rejection  temperature.  This  restriction,  the  Carnot  limitation, 
makes  electrochemical  conversion  a  very  attractive  method.  The  corresponding  set  of 
events  for  a  direct  electrochemical  conversion  is  as  follows:  the  negative  electrode  is 
an  electron  source  while  the  positive  electrode  functions  as  an  electron  sink,  thereby 
causing  a  flow  of  electrons  in  an  external  circuit  without  any  additional  ordering 
process,  and 

-AH^  Uc{i)e-.  (4) 


The  bypass  of  the  Qt  — >  W{pV)  step  enormously  increases  the  conversion  effi¬ 
ciency.  Physically,  the  Qt  term  in  eq.  (3)  represents  the  chaotic  motion  of  molecules 
while  the  term,  Uc{i)e~ ,  is  the  work  done  by  an  orderly  flow  of  electrons  in  an  exter¬ 
nal  circuit.  In  electrochemical  systems,  apart  from  parasitic  chemical  reactions  taking 
place  in  the  cell  interior,  the  chaos  of  the  chemical  reaction  occurs  at  well-defined  po¬ 
sitions  within  the  reaction  space,  i.e.,  at  the  electrode  surface  (Bagotzky  &;  Skundin. 
1980;  Bockris  &  Reddy,  1972). 

The  thermodynamic  efficiency  of  galvanic  cells  is  defined  as  the  ratio  of  electrical 
energy  extracted  to  the  heat  of  reaction.  In  electrochemical  systems  operating  close 
to  the  Uc(0)  ,  the  term  q^nt  =  —T AS  is  negligible.  Therefore,  the  thermodynamic  effi¬ 
ciency  is  nearly  one  and  is  dependent  upon  temperature  through  the  Gibbs-Helmholtz 
relation,  eq.  (5) 


f/o  =  Ut  +  TdUoldT. 
4 


(5) 


Using  this  relationship,  the  thermodynamic  efficiency  may  be  written  in  a  form 
resembling  the  efficiency  of  a  thermal  engine,  viz.  eceii  =  The  temperature- 

dependent  efficiencies  are  compared  in  figure  1.  It  is  seen  that  the  thermodynamic 
efficiency  of  a  LifSOCli  cell  operating  near  the  rest  potential  and  at  the  melting 
point  of  Li  is  ca  0.97.  To  obtain  the  same  efficiency  in  an  ideal  heat  engine,  the 
temperature  of  the  working  fluid  rejected  into  a  300-K  sink,  would  have  to  be  about 
10“*  K.  At  high-discharge  rates,  e.g.,  for  cells  operating  at  80  mA  cm~^,  efficiencies  on 
the  order  of  0.7  were  recorded  experimentally.  Such  a  comparison  might  be  mislead¬ 
ing  inasmuch  as  the  actual  cell  efficiency  is  a  product  of  thermal,  voltaic,  and  faradaic 
efficiencies,  of  which  the  voltaic  and  faradaic  efficiencies  are  temperature-dependent. 
For  practical  reasons,  however,  an  additional  efficiency  can  be  introduced,  namely  an 
engineering  efficiency  that  includes  efficiencies  associated  with  construction  materials, 
reliability  requirements,  cell  design,  parasitic  reactions,  etc.  The  parasitic  reactions 
are  those  that  do  not  contribute  to  the  production  of  electrical  power  and,  at  times, 
may  cause  the  device  to  fail.  The  most  common  examples  of  parasitic  reactions  are 
corrosion,  reactant  and/or  solute  decomposition,  and  film  and  dendrite  formation. 
These  parasitic  reactions  must  be  included  in  the  process  of  battery  optimization  as 
well  as  for  safety  considerations.  In  Li  batteries,  they  arise  because  of  the  highly  re¬ 
active  nature  of  the  cell  components  -  a  prerequisite  associated  with  the  construction 
of  high-energy  and  power-dense  systems. 

Evaluation  of  battery  efficiency  can  be  made  using  an  expression  consisting  of 
two  terms;  the  electrochemical  efficiency  and  the  engineering  efficiency.  Specifically, 
the  first  is  a  dimensionless  number  related  to  the  efficiency  of  all  energy-conversion 
processes,  wherea:  the  second  is  defined  as  the  weight  ratio  of  active  material  to  the 
whole  system,  accessories  included  (Gabano,  1983).  It  should  be  noted,  however, 
there  is  no  general  agreement  as  to  the  procedures  and  test  methods  used  for  system 
evaluation.  For  this  reason  it  is  difficult,  and  often  impossible,  to  compare  results 
reported  by  various  laboratories. 

1.5.  ORGANIZATION  OF  DEVELOPMENT  PROGRAM 

To  illustrate  the  challenges  encountered  during  the  development  phase,  we  con¬ 
sider  the  following  set  of  requirements:  a  battery  voltage  of  300  V  and  power  output 
of  ca  9  kW  for  a  period  of  ca  10  minutes  with  an  activation  time  no  greater  than  30 
seconds.  Moreover,  because  of  space  limitations  and  weight  restrictions,  an  energy 
density  of  560  Wh  dm~^  and  200  Wh  kg'*  w'ould  have  to  be  realized.  The  module 


TEMPERATURE  (°K) 

Figure  1.  Thermodynamic  efficiency  as  a  function  of  temperature,  (a)  Efficiency  of  an 
ideal  heat  engine  with  a  300-K  sink  temperature,  (b)  efficiency  of  a  Lif  SOCI2  cell 
at  indicated  cell  potentials,  V. 


necessary  to  meet  these  requirements  in  a  static  electrolyte  system  must  consists  of  90 
to  100  thin  cells,  of  bipolar  design,  connected  in  series,  each  containing  approximately 
10  to  15  cm^  of  electrolyte.  The  short  activation  time  necessitates  the  installation 
of  more  than  one  filling  port.  The  module  designed  to  satisfy  these  requirements  is 
shown  in  figure  2. 


During  the  development  period,  a  number  of  failures  were  reported.  These  can  be 
divided  into  three  categories,  viz.,  failures  detected  during  module  assembly;  those 
that  occurred  in  the  course  of  module  discharge,  inclusive  of  the  post-discharge  pe¬ 
riod;  and  those  which  led  to  an  unacceptable  performance,  e.g.,  inadequate  battery 
lifetime.  Intercell  short  circuits  were  the  major  problem  during  the  manufacturing 
cycle.  These  undoubtedly  were  due  to  the  contradictory  nature  of  the  system  require¬ 
ments.  The  thin  cell  design  meant  a  separator  thickness  on  the  order  of  5  x  10“^  m 
in  the  compressed  cell.  The  separator  must  be  highly  porous  to  minimize  the  ohmic 
drop  across  the  cell  and  permit  rapid  filling  during  the  activation  process.  In  ad¬ 
dition,  the  separator  must  be  mechanically  strong  to  prevent  distortion  and  tearing 
during  assembly  and  activation  and  must  act  as  a  barrier  to  eliminate  any  migration 
of  carbon  particles.  This  balance  between  porosity  and  strength  leaves  little  margin 
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Figure  2.  An  experimental  9-kW  module  of  thin  cell  design,  reserve  type,  built  by 
Altus  Corp.,  San  Jose,  CA. 


for  error  in  construction.  Figure  3  presents  the  chronology  of  failures  and  corrective 
measures  in  the  test  program.  Most  of  these  failures  were  attributed  to  short  circuits. 
Even  thougli  the  repeated  failures  were  the  consecpience  of  inadequate  quality  control 
procedures,  the  initial  remedy  was  to  increase  cell  thickness,  followed  by  a  change  in 
the  activation  procedure.  The  final  remedy  was  to  revise  the  quality  control  measures 
and  to  improve  manufacturing  techniques. 


Examination  of  the  third  catc'gory.  i.e..  failure  to  meet  the  performance  criteria, 
can  be  traced  to  the  lack  of  a  clear  understanding  of  ti'chnology  issues  and  transition 
activities.  Transition  is  (h'linefl  as  a  passage  from  one'  stage  to  another.  In  the  con¬ 
text  of  battery  development,  the  identifiable  stages  are;  research,  directed  research, 
(hnelopment  process,  etc.  We  define  n’search  as  an  actir  ity  that  seeks  an  extension 
of  knowledge  for  its  own  sak<’:  directed  n'search  is  one  that  examiiu's  observations  for 
the  pur{>ose  of  isolating  tc'chnology  issiu's:  and  the  d('velo|)m<’ut  procf'ss  is  concerned 
with  ix'duction  to  practice.  With  these  <lefinitions  in  mind,  wi-  pr<'sent  an  idealized 
model  for  transition  acti\ities  in  figuix'  la  and  contrast  it  with  one  ap|)li('d  during 
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thr  roust ruct ion  of  a  niodtilo  latod  at  ra  S  l<\\  and  oiiorgy  density  of  ‘200  \\  li 
figure  111.  F,\identlv,  in  tlu'  eoiirst'  of  tlie  t'arly  rU'velopmetit al  activities  it  was  forgot¬ 
ten  that,  for  tlu'  transition  to  occur,  ideas  or  la'snlts  generah'd  within  the  context  of 
activitit's  of  one  stage'  must  satisfy  thi'  constraints  imposed  l\v  t  lu'  lU'xt  stage.  1  lu' 
nature  of  tlu'se'  constraints  may  vary,  l>nl.  in  getu'ial.  tliey  can  Ite  re'duce'd  to  hu¬ 
man  factors,  including  numagi’inent  philosopliv  and  te'chnical  conside'rat  ie»ns  (Szpak 
.k‘  Smith.  !!)!)()).  I'lu'  course  of  this  dc'Ve'lopmenl  |)rogram  jU'ovith's  clear  ('videncc' 
that  the  a[)|)licat  ion  of  solid  enginer-ring  principles  is  le.ss  ('xpe'iisi ve-  and  less  time' 
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consuming  than  intuitive  solutions  and  rapid  fixes. 

2.  CHEMISTRY  AND  ELECTROCHEMISTRY 
OF  THE  LITHIUM/THIONYL  CHLORIDE 

SYSTEM 

2.1.  INTRODUCTORY  COMMENTS 

While  the  energy  density  of  a  system,  i.e.,  the  total  useful  energy  content  is  deter¬ 
mined  by  the  thermodynamics  of  the  galvanic  couple,  its  rate  capabilities  are  deter¬ 
mined  by  electrochemical  kinetics  and  transport  phenomena.  This,  in  turn,  involves 
an  understanding  of  the  partial  reactions  and  processes  that  occur  at  the  individual 
electrodes  and  the  behavior  of  the  electrolyte  during  discharge.  The  thermodynamic 
and  kinetic  data  are  then  used  in  modeling  exercises  to  follow  the  overall  system 
behavior. 

2.2.  LITHIUM/THIONYL  CHLORIDE  OVERALL  CELL 
REACTIONS 

In  determining  the  correct  reaction  stoichiometry,  a  battery  containing  known 
amounts  of  Li  and  SOCI2  is  assembled  and  discharged  at  a  constant  rate.  The 
amount  of  charge  passed  is  monitored  and  when  discharge  is  complete,  the  battery  is 
disassembled  and  examined  for  evidence  of  corrosion  or  dendrite  formation.  Finally, 
the  electrodes  and  electrolyte  are  subjected  to  chemical  analysis  to  determine  the 
products  of  electroreduction. 

Table  1  lists  a  number  of  proposed  overall  cell  reactions  (Szpak  &  Driscoll,  1987). 
The  fact  that  as  many  as  12  reaction  schemes  have  been  proposed  indicates  the 
complexity  of  the  chemistry,  as  well  as  the  experimental  difficulties  encountered  in 
analyzing  this  system.  A  brief  review  of  the  early  literature  indicates  disagreement 
among  investigators.  For  example,  Auborn  et  al.  (1973)  suggested  that  eq.  (II), 
table  1,  adequately  represents  the  cell  reaction,  while  Dey  (1976)  postulated  that 
eq.  (I),  was  the  correct  reaction.  In  1976,  Driscoll,  Holleck,  and  Toland  reported  an¬ 
alytical  data  in  favor  of  eq.  (I),  giving  Li'*',  Cl~  and  S  in  the  expected  stoichiometries. 
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Table  1.  Summary  of  postulated  overall  cell  reactions  for  the  Li/SOCl2  system. 
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Energy  density  calculated  from  total  energy  (-AG,.)  according  to  the  equation 
-.AG^ 

=  ^  -A - r  where  -SGr  is  in  Joules 

3.6(6.94nLi+ tl8.9nsoci2> 

and  nLj  nsoci2  stoichiometric  coefficients  from  column  2 


They  also  identified  SO2  as  a  reaction  product;  however,  the  amount  of  SO2  found  wa^ 
not  consistent  with  eq.  (I).  It  was  later  shown  that  a  side  reaction  in  the  analytical 
procedure  resulted  in  this  discrepency.  In  1979,  Schlaikjer  et  al.  reported  that  the 
reaction  pathway  is  dependent  upon  the  operating  temperature  of  the  cell,  viz.,  at 
temperatures  above  253  K,  eq.  (I),  table  1,  dominates,  while  below  253  K  eq.  (II) 
represents  the  correct  reaction.  Therefore,  on  the  basis  of  the  best  presently  available 
evidence,  eq.  (I),  table  1,  is  considered  to  represent  the  overall  cell  reaction  and  is 
used  as  the  basis  for  the  modeling  of  the  Li/SOCl2  system. 

2.3.  PROPOSED  MECHANISMS  OF  THIONYL  CHLORIDE 
ELECTROREDUCTION  -  AN  EARLY  PERIOD 


The  interest  in  determining  the  mechanism  of  SOCI2  electroreduction  dates  back 
to  the  mid  1970s  and  extends  through  the  1980s.  The  methods  employed  include 
cyclic  voltammetry,  spectroscopic  examination  of  spent  electrolytes,  gas  chromatog¬ 
raphy,  and  electron  spin  resonance  (esr)  spectrometry.  In  all  reported  mechanisms, 
the  key  role  was  assigned  to  the  formation  of  SOCl  radical: 


SOCI2  +  e-  SOCl  +  Cl-  ,  (I) 

which  then  reacts,  via  a  series  of  steps,  to  form  the  products  of  electroreduction.  The 
proposed  reaction  schemes  are  summarized  in  table  2.  Bowden  and  Dey  (1980)  exam¬ 
ined  the  reduction  of  cell  electrolytes  using  cyclic  voltammetry.  To  avoid  difficulties 
associated  with  electrode  passivation  by  the  precipitation  of  LiCl,  they  examined 
voltarnniograms  of  elemental  sulfur,  SOCI2  i-S’O?,  and  S2CI2  dissolved  in  organic  sol¬ 
vents,  e.g.,  A^(C/f3)2C'//0;  CH3CN:  {C  113)280  .  From  the  shape  of  voltammograms, 
augmented  by  ultraviolet- visible  (UV-VIS)  spectra  of  cell  electrolytes  sampled  during 
discharge,  they  concluded  that,  after  the  formation  of  the  SOCl  radical,  the  reduc¬ 
tion  proceeds  via  the  production  of  a  hazardous  intermediate  (S'O)!.  Investigations 
by  Istone  and  Brodd  (1984)  further  supported  the  Bowden  and  Dey  mechanism. 
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Table  2.  Postulated  reaction  mechanisms  for  the  reduction  of  thionyl  chloride 
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Blomgren  et  al.  (1979)  also  took,  as  the  starting  point,  the  formation  of  SOCl 
radical,  followed  by  a  number  of  steps  yielding  SCI2  and  S2CI2  without  the  formation 
of  a  hazardous  intermediate.  In  1985,  Carter  et  al.  (1985)  using  modern  tools,  includ¬ 
ing  esr  spectrometry,  monitored  the  electroreduction  of  LiAlCli! SOCI2  electrolytes. 
Their  gas  chromatographs  showed  an  increase  in  C/2,  SO2,  and  SO2CI2  in  the  dis¬ 
charged  samples.  Similarly,  the  infrared  (IR)  spectra  indicated  formation  of  SO2  in 
the  discharged  solutions  as  well  as  formation  of  S02C12-  An  esr  signal  was  observed 
for  these  solutions,  implying  the  existence  of  a  radical,  presumably  SOCl.  They  pro¬ 
posed  a  rather  complex  mechanism  with  the  formation  of  a  hazardous  intermediate, 
CC/5-,  cf.  table  2. 


It  is  noteworthy  that  the  objective  of  this  research  was  to  find  an  intermediate 
responsible  for  the  catastrophic  events  during  discharge.  For  example,  Salmon  et  al. 
(1982)  reported  formation  of  C/2C  and  CIO2;  however,  later  work  (Carter,  Frank 
&  Szpak,  1984)  failed  to  confirm  these  findings.  At  present,  it  is  believed  that  the 
forced  explosive  venting,  which  sometimes  occurs  during  the  discharge  of  Li/SOCh 
battery  is  due  to  the  development  of  hot  spots  resulting  from  either  internal  shorts 
or  localized  exothermic  reactions  between  the  metallic  lithium  and  SOCI2. 

2.4.  CHEMISTRY  OF  THE  LiCl  -  AlCh  -  SOCh  SYSTEM 


The  proposed  mechanisms  summarized  in  table  2,  ignore  the  fact  that  SOCI2 
forms  adducts  with  A/C/3  (Long  &  Bailey,  1963;  Hecht,  1947;  and  Spandau  h  Brun- 
neck,  1952)  and  is  able  to  enter  the  solvation  sheet  of  ionic  species  either  initially 
present  or  generated  in  the  course  of  the  charge  transfer  process.  If  the  composition 
of  the  electrode/electrolyte  interphase  is  related  to  the  bulk,  eq.  (I)  as  written  does 
not  identify  which  SOCI2  molecules  reside  in  the  reaction  space  and  participate  in  the 
charge  transfer  process.  Furthermore,  Madou  and  Szpak  (1984)  and  Madou,  Szpak, 
and  Smith  (1987)  have  shown,  that  more  than  one  S'OC/2-bearing  species  undergo 
electroreduction.  Consequently,  it  is  essential  to  know  what  species  are  present  in 
solution,  their  structures,  and  equilibria.  Experimental  details  are  given  in  appendix 
A. 
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2.4.1. 


The  Structure  of  Neat  SOCI2 


Thionyl  chloride  is  a  weakly  associated  liquid.  It  has  a  rather  low  latent  heat  of 
vaporization  and  there  is  a  slight  change  in  the  S  —  O  stretching  frequency  due  to 
condensation  (the  S~0  stretching  vibration  occurs  at  1251  cm“^  in  the  gas  phase  and 
1231  cm~^  in  the  liquid  phase).  Further  evidence  of  weak  association  is  the  numerical 
value  of  the  slope  of  the  fluidity  as  a  function  of  specific  volume,  which  indicates 
the  presence  of  small,  interacting  molecular  clusters  (Szpak  &;  Venkatasetty,  1984). 
Examination  of  the  S  —  O  stretching  vibration  of  neat  SOCI2  at  1231  cm“^,  shown 
in  figure  5,  indicates  that  it  is  actually  a  composite  band  composed  of  two  Voigt 
profiles.  When  cooled  to,  e.g.,  250  K  the  band  due  to  the  S  —  O  stretching  vibration 
is  composed  of  three  Voigt  profiles,  (Mosier-Boss  et  al.,  1989a).  These  Voigt  profiles 
have  been  assigned  to  dimers  and  oligomers  of  SOCh-,  which  according  to  molecular 
orbital  (MO)  calculations  have  an  open-chain  spatial  structure,  shown  in  figure  6. 

2.4.2.  Structure  of  the  AI2CIQI SOCh  System 

With  the  addition  of  AhCle,  a  new  set  of  bands  appear  in  the  Raman  spectra, 
as  shown  in  figure  7,  indicating  the  formation  of  adduct(s).  At  low  concentrations  of 
AhCh,  a  new  band  at  1108  cm"'  appears.  This  band  has  been  assigned  to  the  S  —  O 
stretching  vibration  of  the  1:1  adduct  (Long  &  Bailey,  1963).  Since  the  observed 
shift  of  the  S  ~  O  stretching  vibration  occurs  at  lower  frequency  and  the  symmetric 
and  asymmetric  S  ~  Cl  stretching  vibrations  shift  to  higher  frequency,  complexation 
occurs  through  the  oxygen  atom,  ClsO  AICI3.  As  the  AhCle  concentration  in- 
crea.ses,  a  new  S  —  O  stretching  vibrational  band  at  ca  1055  cm"'  appears,  as  shown 
in  figure  7.  This  new  band  is  very  broad  and  is  due  to  a  5  —  O  stretching  vibration  of 
a  second  AICI3  —  SOCI2  species.  The  broadness  of  this  band  as  well  as  the  shift  of  the 
S  —  O  stretching  vibration  to  lower  frequency  indicates  further  electron  withdrawal 
from  the  oxyger  atom  perhaps  from  the  formation  of  a  1:2  complex,  Cl2SO{AlCl3)2, 
as  suggested  by  others  (Long  Bailey,  1963;  Hecht,  1947).  Such  a  conclusion  would 
be  consistent  with  the  existence  of  a  1:2  complex  in  the  solid  state  (Long  Bailey, 
1963);  however,  the  appearance  of  this  band  coincides  with  an  increase  in  the  solution 
conductivity  (Szpak  &  Venkatasetty,  1984).  Molecular  adducts  are  neutral  and,  as 
such,  are  nonconductive.  Therefore,  this  new  band  occurring  at  1055  cm"'  is  due  to 
an  ionic  species  and  not  to  Cl2SO{AlCl3)2-  Formation  of  an  ionic  species  from  the 
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Raman  and  l.R.  Spectroscopy  o/AlCl,-SOa, 
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wavenumber/cm'' 

Figure  5.  Decomposition  of  S  —  0  stretching  band  of  neat  SOCI2  into  Voigt  profiles, 
(a)  273.5  K,  (b)  250.0  K  (Mosier-Boss  et  al.,  1989a). 


CI2SO  AICI3  adduct  can  occur  either  by  a  halide  ion  transfer 


CI2SO  ^  AICI3  ^  SOCI+  +  AICI4  ,  (II) 

or  by  an  internal  exchange  reaction 


2CI2SO  ^  AICI3  -  ChAli^  OSCl2)t  +  AlCl:  .  (Ill) 
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Figure  6.  Molecular  orbital  generated  structures  for:  (a)  {SOCl2)2  open 
chain  dimer,  (b)  CI3AI  OSCI2  1:1  adduct,  (c)  [Cl2Al{<— 

(d)  [Li{*—  03012)2]^  solvated  Lt'*’  ion,  (e)  [Li{*—  O SC  12)2]^[AIC I4]  ion  pair, 
and  (f)  Li{OSCl2,  S02)^  ion  (Mosier-Boss  et  al.,  1991). 


Figure  7.  Evolution  of  Raman  spectra  of  the  AICI3  —  SOCI2  system  as  a  function 
of  AICI3  concentration  as  indicated.  The  S  —  O  stretching  region  occurs  between 
1000-1250  cm~^.  The  symmetric  and  asymmetric  S  —  Cl  stretching  region  occurs 
between  400-600  cm~*  (Mosier-Boss  et  al.,  1989a). 


Both  the  spectroscopic  data  and  MO  calculations  favor  formation  of  an  onium 
ion,  Cl2Al{*—  05012)2  j  via  eq.  (Ill)  (Mosier-Boss  et  al.,  1989o).  The  onium  ion 
has  a  negative  A/f/  =  — 11  kcal  mole"’  compared  to  SOCl'*',  which  has  a  positive 
AHf  =  -t-228.27  kcal  mole"’.  Since  the  S  —  O  bond  of  SOCl'^  has  a  higher  bond 
order  than  that  for  SOCI2,  the  S  —  O  stretching  vibration  of  SOCl'*'  would  occur  at 
a  higher  frequency.  However,  no  high-frequency  bands  were  observed  in  the  S  —  O 
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stretching  region  of  the  Raman  spectra.  In  contrast,  the  S  —  O  stretching  vibration 
of  CUAl^*—  08012)2  is  expected  to  occur  at  a  lower  frequency  than  for  either  SOCI2 
or  the  1:1  adduct.  The  MO-generated  structure  for  the  1:1  adduct  and  onium  ion  are 
shown  in  figure  6. 


2.4.3.  Structure  of  the  LiCl ! AI2CIQI SOCI2  System 

The  effect  of  LiCl  addition  on  the  structural  features  of  the  AICI3  —  SOCI2 
system  is  illustrated  in  figure  8.  With  the  addition  of  LiCl  the  bands  due  to  the 
1:1  adduct,  CI3AI  *—  0SC12-,  and  onium  ion,  Cl2Al[*—  05C/2)J,  disappear.  In  the 
S  —  O  stretching  region  a  new  band  at  1215  cm“'  appears  which  has  been  assigned 
to  solvated  ion.  Again,  the  lowering  of  the  S  —  O  stretching  vibration  of  solvated 
Li^  indicates  coordination  is  through  the  oxygen  atom,  Li{*~  OSCl2)2-  At  least 
two  SOCI2  molecules  are  required  to  solvate  Li'^  as  a  1:1:2  molar  ratio  solution  of 
LiCl  :  A/C/3  •  SOCI2  can  be  prepared.  The  Raman  spectrum  of  this  solution  shows 
no  free  SOCI2  to  be  present  (Mosier-Boss  et  ah,  1989/)). 


Typically,  Li"^  has  a  coordination  number  of  four.  In  the  presence  of  excess 
SOCh',  these  sites  are  likely  filled  wi\,h  SOCU  ’uo'vcules.  However,  in  the  1:1:2 
LiCl  :  AICI3  :  SOCI2  mole  ratio  so'-itic  ,  two  sites  are  occupied  by  SOCI2  and 
the  remaining  two  coordination  sites  can  be  filled  by  ion  pairing  with  AICI4.  While 
SOCI2  is  a  polar  solve  ,  as  s.^own  by  irs  Ir.rg^  d'oole  moment,  it  has  a  low  di¬ 
electric  constant  implying  that  substantial  ion  pairing  should  occur  in  these  elec¬ 
trolytes.  Indeed,  from  conductance  measurements  of  LiAlCU  electrolytes  in  SOCI2, 
Venkatasetty  and  baathoff  (1981)  obtained  a  value  of  1.605  x  10~^  for  the  dissociation 
constant  of  the  ion  pair.  MO  calculations  of  solvated  Li'^ ,  C5C/2)J,  and  the 

ion  pair,  [Li{*—  05C/2)J][A/C/7],  favor  formation  of  the  ion  pair.  Not  only  is  the 
heat  of  firmation.  A///,  more  negative  for  the  ion  pair  (—291  kcal  mole“^  for  the  ion 
pair  and  —56  kcal  mole"'  for  the  solvated  Li'*'  ion),  but  the  Li  —  O  bond  length  for 
the  Li'*'  <—  OSCI2  interaction  decreases  and  the  positive  charge  gets  distributed  into 
the  SOCI2  moieties  of  the  ion  pair.  The  MO-generated  structures  for  Li{*~  OSCl2)t 
and  ion  pair  are  shown  in  figure  6. 


2.4.4.  A  Brief  Summary  of  the  SOCI2  Chemistry 


The  structures  and  stiuctural  changes  of  these  electrolytes  can  be  summarized 
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RAMAN  FREQUENCY  (cm'^ 


Figure  8.  Raman  spectra  of  the  LiCl  —  AICI3  —  SOCI2  system.  3.0  M  AICI3 
in  SOCh.  (a)  0  M  LiCl,  (b)  1.0  M  LiCl,  (t)  2.0  M  LiCl,  and  (d)  3.0  M  LiCl 
(Mosier-Boss  et  al.,  19896). 


as  follows.  In  the  liquid  state,  SOCI2  forms  open-chain  dimers  and  oligomers.  The 
dissolution  of  AI2CIQ  occurs  dissociatively  with  the  formation  of  adducts  which,  at 
higher  Al2Cle  concentrations,  yield  onium  ions  and  AlCl^.  The  addition  of  LiCl 
causes  the  destruction  of  CI3AI  OSCI2  and  the  onium  ion  with  the  formation  of 
solvated  Li'^ .  At  higher  LiAlCU  concentrations,  ion  pairing  between  the  solvated 
Li+  ion  and  AlCl^  occurs.  The  spectral  data,  results  of  the  MO  calculations,  and 
the  viscosity  and  conductivity  measurements  are  in  agreement.  These  results  indicate 
that  the  S  —  O  bond  strength  of  the  free  and  complexed  SOCI2  diminishes  as  follows: 

CI2SO  >  (CI2SOU  >  Li{*-  OSCl2)t  > 

ChAl  ^  OSCI2  >  ChAli^  OSCl2)t;  n>2. 
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2.5. 


ELECTROCHEMICAL  ASPECTS  OF  SOCh  REDUCTION 


The  reaction  path  of  electroreduction  of  SOCI2  electrolytes  is  examined  by  three 
techniques:  in  situ  reflectance  infrared  spectroscopy,  cyclic  voltammetry,  and  both 
galvanostatic  and  potentiostatic  pulsing.  Experimental  details  are  given  in  appendix 
A.  What  follows  is  a  brief  discussion  of  the  techniques  chosen  to  probe  these  electro¬ 
chemical  reactions. 


2.5.1.  In  Situ  Reflectance  Infrared  Spectroscopy 


In  situ  reflectance  infrared  spectroscopy  is  a  useful  technique  to  identify  the  species 
present  in  the  electrode-electrolyte  interphase  region  (Ashley  Sz  Pons,  1988).  This 
technique  provides  information  regarding  surface  coverage,  the  orientation  of  adsor¬ 
bates,  as  well  as  structural  information  on  the  intermediate  and  product  species 
formed  upon  electroreduction.  By  varying  the  electrode  potential,  it  is  possible  to 
determine,  from  changes  in  the  spectra,  the  species  undergoing  electroreduction  a.s 
well  as  their  orientation  on  the  electrode  surface. 

2.5.2.  Cyclic  Voltammetry 


Cyclic  voltammetry  is  the  experimental  technique  of  choice  to  assess  the  complex¬ 
ity  of  an  electrochemical  reaction  and  elucidate,  at  least  qualitatively,  its  mechanism. 
The  shape  of  the  current/voltage,  alternatively,  current/time  curves,  reflects  the  spe¬ 
cific  set  of  events  associated  with  the  charge  transfer  reaction  taking  place  under 
selected  experimental  conditions.  The  current/time  relationship,  i(0,t),  for  a  linear 
potential  scan,  U{t)  =  U°  +  |ult,  for  the  N  electroactive  species  and  M  adsorption 
processes,  is  given  by 


j(0,t)  =  F 


N 


M 


d0„ 


Tn=l 


(^>) 


where  denotes  the  maximum  surface  concentration  of  the  m  —  th  species  and  6 
is  the  surface  coverage.  The  right-hand  side  of  eq.  (6)  contains  contributions  due  to 
the  diffusional  flux  via  the  sum  of  concentration  gradients  of  electroactive  species  and 
the  kinetics  of  adsorption  expressed  in  terms  of  the  time  rate  of  change  in  the  surface 
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coverage  summed  over  the  number  of  adsorption  processes.  Which  set  of  processes 
dominates  the  charge  transfer  reaction  depends,  among  other  factors,  on  the  rate  of 
linear  scan.  In  general,  slow  scan  rates  relate  to  the  diffusional  processes  while  fast 
scan  rates  expose  the  adsorption  characteristics.  The  concentration(s)  gradient  of  the 
electroactive  species  is  obtained  from  the  mass  balance  equation,  eq.  (7) 

|  =  (7) 

which  is  solved  subject  to  boundary  conditions  reflecting  the  experimental  constraints. 
The  function,  fni(,  A)  ,  is  formulated  for  a  specific  set  of  events  occurring  within  the 
interphase  in  the  course  of  charge  transfer.  In  particular,  ^  refers  to  thermodynamic 
aspects,  e.g.,  equilibria  between  species  populating  the  interphase,  while  A  contains 
information  on  the  rate  of  relevant  processes,  including  forcing  perturbation,  e.g., 
potential  scan  rate. 

2.5.3.  Pulsing  Techniques 


Additional  information  can  be  extracted  by  changing  the  nature  of  the  forcing 
function,  /n($,  A).  For  a  class  of  electrode  processes  where  adsorbed  species  dominate, 
the  electrode  potential  can  be  written  in  a  general  form, 

u  =  (8) 


and  its  time  dependence  can  be  obtained  by  formal  differentiation  which,  for  the 
galvanostatic  case  yields 

dt  “  V  dt  ’ 


Galvanostatic  and  potentiostatic  pulsing  techniques  are  often  used  to  probe  the 
kinetics  of  mass  transport,  charge  transfer,  coupled  chemical  reaction  rates  and  ad¬ 
sorption  (Murray,  1986;  Bard  &  Faulkner,  1980).  In  the  galvanostatic  pulsing  exper¬ 
iment,  the  constant  current  applied  to  the  electrode  causes  the  electroactive  species, 
ox,  to  be  reduced  to  product,  red.  This  process  is  essentially  a  titration  of  the  ox 
species  near  the  electrode  by  the  constant  flux  of  electrons.  The  potential  of  the 
electrode  moves  to  values  characteristic  of  the  couple  and  varies  with  time  as  the 
Cox/Cred  ratio  changes  at  the  electrode  surface  resulting  in  a  titration  curve.  In  the 
potentiostatic  pulsing  experiment,  a  potential  step  is  applied  to  the  working  electrode 
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and  the  resulting  current/time  response  is  measured.  The  magnitude  of  the  current 
at  any  given  time  and  the  shape  of  the  current/time  response  are  also  dependent 
upon  the  mass  balance  equation  through  the  forcing  function,  /n(^,  A),  eq.  (7). 


2.6.  ELECTROREDUCTION  OF  SOCI2  ELECTROLYTES  - 
REACTION  PATH 


As  expected,  a  given  electrolyte  composition  may  be  ideal  for  one  battery  ap¬ 
plication  but  not  another.  Electrolyte  composition  strongly  influences  the  reaction 
path  of  SOCI2  electroreduction.  In  practice,  the  high-energy  dense,  reserve  batteries 
usually  employ  acidic  electrolytes  whereas  the  low-rate  cells  use  neutral  electrolytes. 
It  was,  therefore,  of  interest  to  us  to  determine  the  similarities  and  differences  in  the 
electroreduction  pathways  of  acidic  and  neutral  electrolytes. 

2.6.1.  Concept  of  an  Electrode— Electrolyte  Interphase 


An  interphase  region  is  formed  whenever  an  electrode  is  in  contact  with  an  elec¬ 
trolyte.  In  the  simplest  case,  the  interphase  region  takes  the  form  of  the  electrical 
double  layer.  In  the  more  complex  cases  and,  in  particular,  during  the  charge  transfer 
reaction,  it  consists  of  layers,  each  associated  with  a  participating  elementary  process. 
In  this  representation,  the  interphase  region  is  an  open  system  in  which  a  number 
of  consecutive  processes  takes  place,  of  which  the  slowest  one  determines  the  rate. 
These  processes  include  transport  of  the  reactants  from  the  bulk  to  the  electrode 
surface  by  diffusion,  adsorption  on  the  electrode  surface,  charge  transfer,  desorption 
of  the  reaction  products,  followed  by  transport  of  the  reaction  products  away  from 
the  electrode  surface.  In  a  discharging  battery,  these  same  processes  occur;  however, 
in  a  battery  the  electrons  ultimately  flow  into  an  external  circuit  where  the  electrical 
work  is  delivered. 

2.6.2.  Interphase  at  Rest 


Figure  9  shows  the  S  —  O  stretching  region  of  the  AICI3  —  SOCI2  system  re¬ 
flected  off  gold  at  the  UdO).  Varying  the  path  length  between  the  Au  electrode 
and  the  infrared  transparent  window  clearly  indicates  an  enrichment  of  both  the  1;1 
adduct  and  onium  ion  near  the  electrode.  Deconvolution  of  the  composite  band  in  the 
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950-1150  cm“^  infrared  spectral  region  into  Voigt  profiles,  shown  in  figure  10,  pro¬ 
vides  additional  information  on  the  structure  of  the  interphase  at  the  rest  potential. 
The  band  in  the  transmittance  spectrum  consists  of  two  Voigt  profiles  at  1116.9  and 
1065.5  cm~^,  which  are  assigned  to  the  adduct  and  onium  ion,  respectively.  However, 
the  reflectance  spectra  of  Au  and  Pi  show  two  additional  Voigt  profiles  at  1132  and 
1102  cm"^.  These  two  peaks,  which  disappear  upon  cathodic  polarization,  were  ten¬ 
tatively  assigned  earlier  as  either  (i)  the  symmetric  and  asymmetric  S  —  O  stretching 
vibrations  of  adsorbed  onium  ion  or  (ii)  the  S  —  O  stretching  vibrations  of  adsorbed 
onium  ion  and  1:1  adduct,  respectively  (Mosier-Boss  et  al.,  19896).  The  disappear¬ 
ance  of  the  1132  and  1102  cm~^  peaks  upon  cathodic  polarization  appears  to  be  cou¬ 
pled,  suggesting  they  are  derived  from  the  same  species.  Therefore,  the  most  probable 
assignment  of  these  two  peaks  is  to  adsorbed  onium  ion.  Regardless,  since  these  peaks 
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Figure  9.  IR  reflectance  spectra  at  rest  potential  at  path  lengths  of  5  and  10  fim, 
respectively.  Electrolyte  4.0  M  AICI3  in  SOCI2,  electrode  material  Au  (Mosier- 
Boss  et  al.,  1989c). 
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are  shifted  to  higher  frequencies,  it  follows  that  adsorption  occurs  through  the  sulfur 
atom.  This  conclusion  is  supported  by  MO  calculations,  which  indicate  that  the 
positive  charge  rests  on  the  sulfur  atom  (Mosier-Boss  et  al.,  1989a).  Therefore,  at 
rest  potential  the  interphase  of  the  MefAlCh  —  SOCI2  system,  illustrated  in  figure 
11,  consists  of  an  electrode  in  contact  with  an  adsorption  layer,  which  is  in  contact 
with  an  enrichment  layer  of  1:1  adduct,  Ai,  and  onium  ion  ,  which  is  in  contact 
with  the  bulk  electrolyte. 
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Figure  10.  The  950-1150  cm~'  IR  spectral  region  decomposed  into  Voigt  profiles, 
(a)  Transmittance  spectrum,  (b)  reflectance  spectrum  from  Au,  (c)  reflectance  spec¬ 
trum  from  Pt.  Both  reflectance  spectra  are  at  the  rest  potential.  A,C,  S  —  O 
vibrations  due  to  onium  ion  and  complex  in  the  bulk  electrolyte;  B,D,  asymmetric  and 
symmetric  S  —  O  stretching  vibrations  of  adsorbed  onium  ion,  respectively,  or  S  —  O 
stretching  vibrations  of  adsorbed  onium  ion  and  1:1  adduct,  respectively.  Electrolyte 
4.0  M  /I/C/3  in  SOCI2  (Mosier-Boss  et  al.,  1989c). 
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Figure  11.  Structure  of  the  electrode/electrolyte  interphase  at  the  rest  potential: 
(b)  -  bulk,  (i)  -  interphase,  (a)  -  adsorption  layer  (Mosier-Boss  et  al.,  1989c). 


2.6.3.  Cathodically  Polarized  Interphase 


The  effect  of  applied  overpotential  on  the  S  —  O  stretching  region  of  Au/4.0  M 
AlClz  in  SOCI2  is  shown  in  figure  12.  With  the  passage  of  current,  new  peaks  appear, 
viz,  at  1331  and  ca  1150  cm~^  due  to  the  formation  of  SO2  and  at 
1170  cm"^  and  1190  cm~^  assigned  to  intermediate  species  Pi  and  P2  with  the  peak  at 
1170  cm"*  appearing  first.  Concurrently,  there  is  a  loss  of  the  peaks  due  to  the  Ai  and 
Aj  species.  These  results  indicate  that  the  electroactive  species  is  not  neat  SOCI2, 
Aq.  Instead  the  species  undergoing  electroreduction  is  either  the  1:1  adduct,  Ai,  or 
the  onium  ion,  At,  or  both.  It  is  not  clear  whether  the  SOCI2  moieties  of  the  1:1 
adduct  and  onium  ion  are  simultaneously  reduced  or  if  only  one  species  participates 
in  the  reduction  process.  This  ambiguity  exists  because  the  species  are  coupled  to  one 
another  via  the  chemical  equilibrium,  eq.  (Ill),  and  products  of  the  electroreduction 
of  one  species  may  shift  the  equilibrium  to  regenerate  the  electroreducible  species. 
In  these  potential-dependent  spectra,  there  is  no  evidence  of  SO,  {SO)x,  SOCl  or 
SO2CI2  intermediate  species  a.s  proposed  by  others,  vide  supra. 


2.6.4.  Relaxation  Techniques 


A  useful  technique  to  explore  complex  electrochemical  processes  is  cyclic  voltam¬ 
metry,  since  the  shape  of  the  current/ voltage  curves  depends  upon  the  thermodynam¬ 
ics  of  the  system  and  the  kinetics  of  the  charge  transfer  and  associated  transport  pro¬ 
cesses,  eqs.  (6)  and  (7).  Cyclic  voltammograms  obtained  for  the  PtjAlCls  —  SOCI2 
systems  are  shown  in  figure  13.  For  3.0  M  AIC/3  in  SOCI2,  two  electron  transfers 
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Figure  12.  The  900-1500  cni~'  IR  spectral  region  of  cathodically  polarized  Au/A.Q  M 
AlClz  in  SOClz-  Overpotentials  vs  Ag/AgCl:  (a)  0.0  V,  (b)  —0,5  V,  (c)  —1.0  V, 
(d)  —1.5  V,  (e)  —2.2  V  (Mosier-Boss  et  al.,  1991). 


are  observed  of  which  the  first,  to  form  intermediate  species  Pj,  is  irreversible  and 
the  second  electron  transfer,  to  form  P2,  is  quasi-reversible.  A  current  plateau  is  as¬ 
sociated  with  the  first  electron  transfer  suggesting  a  coupled  chemical  reaction  either 
preceding  or  parallel  with  the  charge  transfer  (Nicholson  h  Shain,  1964).  Indeed, 
an  equilibrium  does  exist  between  onium  ion  and  1:1  adduct  and,  from  the  in  situ 
infrared  spectroscopy,  a  loss  of  both  these  species  is  observed  with  increasing  over¬ 
potential.  As  the  scan  rate  is  decreased  (figure  13a)  or  the  concentration  of  AlClz 
is  increased  (figure  13b),  an  autocatalytic  effect,  manifested  by  a  cross  over  behavior 
in  the  cyclic  voltammograms  (Nicholson  &  Shain,  1964),  is  observed.  Increasing  the 
AlClz  concentration  increases  the  concentration  of  onium  ions.  Moreover,  since  the 
cross-over  point  occurs  at  the  juncture  between  the  current  plateau  and  the  cathodic 
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Figure  13.  Effect  of  scan  rate  and  concentration  on  voltammogram  shape  for  the 
Pt/AlCI;^  —  SOCI2  system,  (a)  3.0  M  AICI3,  solid  line  -  10  mV  s“^,  dashed 
line  -  5  mV  s~*,  (b)  4.0  M  AICI3,  scan  rate  10  mV  s  *  (Mosier-Boss  et  al.,  1991). 
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peak  due  to  the  second  electron  transfer,  the  species  being  regenerated  on  the  elec¬ 
trode  surface,  via  eq.  (IV),  is  P\. 


P2  -t-  At  ^  2Pi  .  (IV) 

These  results  indicate  that  the  electroactive  species  is  the  onium  ion  and  that 
intermediate  species  Pi  and  P2  are  actually  onium  ion  that  has  accepted  one  and 
two  electrons,  respectively.  The  MO  calculations  show  that  neutral  onium  ion  and 
negatively  charged  onium  ion  are  more  stable  than  positively  charged  onium  ion 
(Mosier-Boss  et  al,  19896). 


At  +  e-  ^  A2,  AHf  =  -310  kcal  mole"^  (V) 

A2  -f  e~  — >  A^iAHf  =  —347  kcal  mole"*.  (VI) 

Moreover,  since  AH/  for  At  is  only  —11  kcal  mole"*,  one  would  expect  the  first 
electron  transfer  to  be  irreversible  while  the  second  would  be  reversible. 


Galvanostatic  ai.d  ^^otentiostatic  pulsing  experiments  provide  additional  infor¬ 
mation  on  the  nt  xiiber  and  charge  of  adsorbed  species.  Figure  14a  shows  the  poten¬ 
tial/time  behavior  across  a  charging  and  relaxing  interphase  for  the  Pi [ AICIZ—SOCI2 
system.  Three  exponential  decays  are  observed  in  the  relaxation  curve.  In  general, 
the  number  of  exponential  decays  present  in  the  potential  relaxation  curve  reflects  the 
num*jer  of  adsorbed  species  (Conway  &  Dzieciuch,  1963).  Therefore,  there  are  three 
adsorbed  species  present  on  the  Pt  electrode.  The  current/time  behavior  during  and 
after  a  potentiostatic  pulse  is  applied  is  shown  in  figure  14b.  Only  two  exponential 
decays  are  observed  in  the  relaxation  curve.  Since  only  charged  species  are  conduc¬ 
tive,  two  of  the  absorbed  species,  in  particular  At  and  P21  are  charged.  Again,  this 
supports  the  notion  that  Pi  and  P2  are  actuallj  the  neutral  and  negatively  charged 
onium  ion,  respectively. 

2.6.5.  Summary  of  Events;  Acidic  Electrolyte 


The  sequence  of  events  occurring  during  SOCI2  electroreduction  in  acidic  elec¬ 
trolytes  is  summarized  in  figure  15.  In  figure  15,  ^  reflects  the  thermodynamics  and 
A  the  kinetics  of  the  system.  The  rate  constants  for  adsorption,  autocatalysis  and 
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Figure  14.  (a)  Potential/time  behavior  and  (b)  current/time  behavior  across  a  charging 
and  relaxing  interphase.  System:  Pt j SOCI2—  2.0  M  AICI3  (Mosier-Boss  et  al., 
1991). 


desorption  are  denoted  by  ka,  k^"^^  and  k^,  respectively,  whereas  the  rate  constants 
for  the  first  and  second  electron  transfers  are  given  by  k^^^  and  k^^\  Surface  cov¬ 
erages  of  Af,  Pi  and  P2  are  defined  by  Oi,  62  and  63,  respectively.  According  to 
this  scheme,  the  intermediate  species.  Pi  and  P2,  desorb  from  the  electrode  surface 
into  the  reaction  layer  where  they  react  with  other  species  to  form  the  products  of 
electroreduction  -  Cl~,  SO2,  and  S.  This  reaction  layer  is  sandwiched  between  the 
enrichment  and  diffusion  layers.  It  is  interesting  to  note  that  in  the  charging  poten¬ 
tial/time  curve  shown  in  figure  14a,  the  potential  does  not  rapidly  fall  off  to  more 
negative  potentials  indicating  that  the  electroactive  species  is  continuously  being  re¬ 
generated  in  the  reaction  layer. 


In  situ  infrared  spectroscopy  of  acidic  electrolytes  from  Au,  Pt,  n  —  Si  and 
glassy  C  surfaces  at  rest  potential  and  cathodically  polarized  are  similar.  From 
this  it  can  be  concluded  that  the  reduction  path  of  acidic  electrolytes,  shown  in  figure 
15,  is  the  same  for  each  surface  (Mosier-Boss  et  al.,  1990).  However,  the  shapes  of  the 
cyclic  voltammograms  and  the  potential/time  and  current/time  curves  upon  charging 
and  relaxation  are  greatly  dependent  upon  the  electrode  material  (Mosier-Boss  et  al., 
1990).  This  is  not  surprising  since  the  shapes  of  these  curves  ultimately  depend  upon 
the  mass  balance  equation,  eq.  (7),  through  the  forcing  function  /n(^,  A).  Therefore, 
the  values  of  the  rate  constants  for  the  processes  summarized  in  figure  15  should  be 
strongly  dependent  upon  the  electrode  material. 
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Figure  15.  Structure  of  the  electrode/electrolyte  interphase  and  summary  of  partici¬ 
pating  processes  (Mosier-Boss  et  al.,  1991). 


2.6.6.  Reaction  Path  in  the  Presence  of  LiCl 


2.6.6. 1.  Interphase  at  Rest  and  Cathodically  Polarized.  Figure  16  shows 
the  in  situ  reflectance  infrared  spectra  of  the  AU/SOCI2  —  AICI3  —  LiCl  interphase 
at  the  rest  potential  and  cathodically  polarized.  At  rest  potential,  a  new  band  due 
to  08012)2  is  observed.  At  cathodic  potentials  less  than  0.9  V,  this  band 

increases  in  size,  indicating  the  accumulation  of  Li{*—  OSCl2)2  in  fhe  interphase  re¬ 
gion.  Formation  of  SO2  is  observed  as  evidenced  by  the  bands  at  1331  and  1150  cm“’. 
the  asymmetric  and  symmetric  S  —  O  stretches  of  SO2  respectively.  As  the  amount 
of  SO2  in  the  interphase  region  increa.ses,  a  new  band  at  1157  cm~'  appears,  which 
has  been  assigned  to  the  symmetric  S  —  O  stretch  of  502  in  the  solvation  shell  of 
Li'^,  i.e.  Li{OSCl2,  802)'^  (Lanson  h  Descroix,  1984).  The  MO-generated  structure 
for  this  species  (AH/  —  —36.4  kcal  mole“*)  is  shown  in  figure  6.  Concurrently, 

there  is  a  loss  of  the  bands  due  to  the  1:1  AICI3  —  8OCI2  adduct  and  onium  ion. 
Therefore,  at  cathodic  overpotentials  le.ss  than  0.9  V,  the  formation  of  5O2  in  elec¬ 
trolytes  containing  both  LiAlCU  and  AICI3  occurs  through  the  electroreduction  of 
onium  ion.  Unfortunately,  the  band  due  to  the  8  —  0  stretch  of  08012)2 
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obscures  the  bands  at  1170  and  1190  cm“*  assigned  to  intermediate  species  P\  and 
P2  (figure  12).  Consequently,  it  is  not  possible  to  assess  the  degree  of  interaction, 
if  any,  between  Pi  and/or  P2  and  Li{*—  OSCl2)2  Li{OSCl2,  S02)'^  present  in 
the  uiterphase  region.  At  cathodic  overpotentials  greater  than  0.9  V,  a  loss  of  the 
band  due  to  Li{*—  O SC  12)2  is  observed  with  further  formation  of  S02-  No  bands 
due  to  Li{S02)t  (Lanson  &  Descroix,  1984)  were  observed  in  the  spectra.  Therefore, 
in  practical  electrolytes,  onium  ion  undergoes  electroreduction  first  followed  by  the 
solvated  lithium  ion.  It  can  be  concluded  that  the  order  of  the  charge  transfer  is 
inversely  related  to  the  strength  of  the  S  —  O  bond.  Similar  changes  in  the  in  situ 
IR  reflectance  spectroscopy  for  the  Pt/AlClz  —  LiCl  —  SOCI2  system  are  observed 
implying  that  the  reaction  paths  are  the  same  for  Au  and  Pt  electrodes. 


Figure  16.  The  900-1500  cm  '  IR  spectral  region  of  cathodically  polarized 
All/  4.0  M  AICI3  and  0.5  M  LiCl  in  SOClz-  Overpotentials  vs  AgJAgCl'. 
(a)  0.0  V,  (b)  -0.6  V,  (c)  -0.9  V,  (d)  -1.5  V  (Mosier-Boss  et  al.,  1991). 
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2. 6. 6. 2.  Relaxation  Techniques.  The  addition  of  a  small  amount  of  LiCl  to 
the  Pt/AlCls  —  SOCI2  system  was  found  to  greatly  modify  the  reaction  path  of 
AICI3  —  SOCI2  electroreduction  as  shown  in  figure  17.  With  the  addition  of  0.1 
M  LiCl  to  the  Pt/  3.0  M  AICI3  in  SOCI2  system,  a  new  cathodic  peak  due  to 
Li(<—  05012)2  electroreduction  is  observed  in  the  cyclic  voltammogram  shown  in 
figure  17a.  By  comparing  the  cyclic  voltammograms  for  that  of  the  Pt / AICI3  —  SOCI2 
system,  figure  13a,  with  that  for  the  PljAlCl^  —  LiCl  —  SOCI2  system,  figure  16a,  it 
can  be  seen  that  the  reoxidation  of  P2  is  suppressed  in  the  presence  of  LiCl  as  is  the  au- 
tocatalytic  reaction,  eq.  (IV).  The  presence  of  LiCl  also  dramatically  suppresses  the 
catalytic  reactions  occurring  in  the  reaction  layer  as  shown  by  the  charging  and  relax¬ 
ation  carves  obtained  from  the  galvanostatic  and  potentiostatic  pulsing  experiments, 
figures  18a  and  18b,  respectively.  In  the  charging  potential/time  curve  shown  in 
figure  18a,  the  potential  is  observed  to  rapidly  fall  off  to  more  negative  values,  indi¬ 
cating  the  flux  of  the  electroactive  species  cannot  keep  up  with  the  electrons  being 
forced  across  the  electrode/electrolyte  interface.  This  contrasts  the  potential/time 
response  observed  in  the  Pt  f  AICI3  —  SOCI2  system,  figure  14a.  Furthermore,  it  takes 
significantly  longer  for  the  interphase  of  the  Pt/AlCl3—LiCl—SOCl2  system  to  return 
to  rest  conditions  after  charging  as  shown  in  figure  18a.  These  results  imply  that  the 
products  of  L?(<—  OSCl2)2  electroreduction,  i.e.,  LiCl  and  Li{OSCl2,  502)'^ ,  act  as 
a  barrier  to  the  diffusion  of  onium  ion  and  Li{i—  OSCl2)2  from  the  bulk  electrolj^te. 
It  is  well  known  that  LiCl  is  not  soluble  in  SOCI2  and  will  precipitate  on  the  electrode 
surface  forming  a  passivating  layer.  However,  in  the  presence  of  CI3AI  <—  OSCI2  and 
onium  ion,  the  LiCl  will  be  scrubbed  off  the  electrode  surface  to  form  AICI4  and 
Li{*—  05012)2  ■  This  may,  in  part,  explain  why  it  takes  longer  for  the  interphase  to 
return  to  rest  conditions  after  charging. 


2. 6. 6. 3.  Summary  of  Events;  Practical  Electrolytes.  Figure  17b  shows  a 
cyclic  voltammogram  obtained  for  a  neutral  electrolyte  off  glassy  C .  Similar  cyclic 
voltammograms  have  been  obtained  for  Pf,  Au  and  n  —  5i  electrodes.  These  results 
indicate  that  the  electroreduction  of  Li{*—  050^)2  is  irreversible.  This  is  supported 
by  MO  calculations.  The  AHj  of  Li{*—  05Cl2)2  is  —56  kcal  mole~^  whereas  for 
the  neutral  species  A///  is  —211  kcal  mole~*.  Therefore,  one  would  expect  the  first 
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Figure  17.  Cyclic  voltammograms  of  the  AICI3  —  LiCl  —  SOCI2  system,  (a) 
System:  Pt/  3.0  M  AlClz  and  0.1  M  LiCl  in  SOCW,  u  =  10  mV  s~^, 
(b)  Glassy  C j  3.0  M  AICI3  and  3.0  M  LiCl  in  SOCW,  v  =  10  mV  s~*  (Mosier- 
Boss  et  al.,  1991). 


Figure  18.  (a)  Potential/time  behavior  and  (b)  current/time  behavior  across  a  charging 
and  relaxing  interphase.  System;  PtfSOCl2—^  0  M  AICI3—Q.I  M  LiCl  (Mosier- 
Boss  et  al.,  1991). 
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electron  transfer,  shown  in  eq.  (VII),  to  be  irreversible. 

Li{^  0SCl2)t  +  e-  ^  2L^(^  OSCh)^  ■  (VII) 

The  MO  calculations  for  the  second  electron  transfer  would  not  converge  indicating 
this  species  is  extremely  unstable.  Apparently,  as  soon  as  the  electron  is  transferred 
to  the  neutral  the  species  falls  apart  into  the  products  of  electroreduction  on  the  elec¬ 
trode.  This  may  contribute  to  the  passivation  observed  in  Li/SOCl2  cells.  Further¬ 
more,  there  is  no  indication  in  the  cyclic  voltammograms  of  autocatalytic  activities 
on  or  near  the  electrode  surface  as  was  observed  for  the  Me/AlCls  —  SOCI2  system. 


3.  MODELING  EXERCISES 

3.1.  INTRODUCTORY  COMMENTS 


The  operation  of  a  battery,  module,  cell,  or  even  a  selected  functional  element 
is  a  dynamic  event  in  which  the  participating  processes  occur  within  a  well-defined 
reaction  space.  This,  in  turn,  implies  the  existence  of  constraints  on  the  individual 
processes.  These  processes  are  interactive  and  can  be  described  by  a  set  of  precise 
mathematical  statements.  In  practice,  the  difficulty  lies  in  formulating  an  acceptable 
model  that  reflects  physical  reality  and,  yet,  at  least  in  the  preliminary  stage,  yields 
solutions  in  a  closed  form.  This  requirement  is  met  only  in  a  few  isolated  cases, 
primarily  because  of  the  nonlinearity  of  the  function  relating  current  to  overpotential. 


Modeling  can  describe  a  large  number  of  operational  situations.  For  the  purpose 
of  illustration,  we  have  selected  several  examples  comparing  computed  results  with 
experimental  data.  Those  presented  here  were  selected  because  of  their  concerns  with 
the  cell  performance  and  safety  of  operation  that  were  pertinent  to  the  LifSOCli  bat- 
tc  y  development  effort  of  the  1980s.  They  are  also  limited  to  the  battery  designed 
for  high-discharge  rates  and  operating  under  conditions  that  tend  to  maximize  the 
extractable  energy  content  at  high-power  output.  The  conditions  addressed  are  elec¬ 
trochemical  processes;  thermal  management,  including  catastrophic  thermal  runaway; 
and  intercell  currents  and  their  respective  effects  on  the  performance  and  safety  of  a 
module  and  assembly  of  modules. 
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3.2.  Li/SOCh  CELL  -  SCHEMATIC 


Realistic  modeling  requires  information  on  the  specification  of  the  mode  of  oper¬ 
ation  and  the  description  of  the  physical  constraints  under  which  the  battery  oper¬ 
ates.  To  illustrate,  section  of  a  Li/SOCh  battery  of  bipolar  construction  is  shown 
in  figure  19.  Each  cell  contains  four  functional  elements:  the  negative  electrode, 
the  positive  electrode,  the  electrolyte,  and  a  porous  separator.  The  separator  serves 
as  a  reservoir  for  the  electrolyte  and  provides  an  ionically  conductive  path.  A  thin 
metal  foil  functions  simultaneously  as  a  chemical  barrier  and  an  electronic  conductor 
between  cells.  Processes  that  occur  in  an  operating  cell  are  diffusion  of  the  participat¬ 
ing  species  toward  the  electrode  surface,  adsorption  on  the  electrode  surface,  charge 
transfer,  desorption  of  reaction  products  and  diffusion  of  the  products  away  from  the 
electrode  surface.  The  flow  of  electrons  in  an  external  circuit  completes  the  process. 


3.3.  CELL/BATTERY  MODELING  -  ELECTROCHEMICAL 
PROCESSES 


The  first  comprehensive  modeling  of  the  dynamics  of  the  discharging  LijSOCl^ 
has  been  presented  by  Tsaur  and  Pollard  (1984a  and  19846).  In  their  model,  the 
mathematical  statements  describing  the  material  balances,  flux  relations  in  the  elec¬ 
trolyte  phase,  as  well  as  the  Faraday’s  and  Ohm’s  laws,  are  unambiguous.  However, 
the  formulation  of  the  governing  equations  involving  electrode  polarization  and  the 
active  specific  surface  in  the  positive  electrode  require  more  care.  The  exponential 
form  of  the  kinetic  relation  usually  employed  is  the  Butler- Volmer  equation 

■sj  =  -sio  [exp(aaF’»?/i?r)  -  ex^{-acFr}l RT)]  (10) 

where  s  is  the  specific  surface  area  and  jo  is  the  exchange  current  density.  This  form 
of  the  kinetic  relation  applies  only  within  a  rather  narrow  range  of  overpotentials  and 
selected  group  of  electrochemical  systems.  Other  forms  of  kinetic  expressions  better 
represent  the  current/potential  relation  of  the  reduction  of  SOCh-  Nonetheless, 
eq.  (10)  has  been  used  because  its  exact  form  does  not  substantially  alter  the  predicted 
trends;  however,  the  numerical  values  of  the  parameters  s  and  a's  are  essential  to 
evaluate  the  cell  lifetime. 


Because  of  the  profound  effect  that  the  specific  surface  area  exerts  on  the  cell 
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Figure  19.  Cross  section  of  a  Lij SOCli  cell. 


lifetime  and  because  it  can,  to  a  certain  degree,  be  manipulated  experimentally,  it  is 
necessary  to  assign  a  functional  form  relating  it  to  the  amount  of  charge  transferred 
locally.  Initially,  an  empirical  expression,  which  assumes  instant  precipitation  of  LiCl 
crystallites  resulting  in  a  permanent  passivation  of  the  electrode  surface,  was  used 


(11) 


where  e  is  the  void  fraction,  v  is  the  volume  fraction  of  the  precipitate  and  the  expo¬ 
nent  p{l,  t)  is  related,  in  general  terms,  to  the  morphology  of  the  LiCl  cystallites.  The 
assumption  of  an  instantaneous  precipitation  of  LiCl  is,  however,  doubtful  although 
good  agreement  has  been  observed,  in  particular  concerning  the  optimum  electrolyte 
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composition  and  the  thickness  of  the  positive  electrode.  Later,  Tsaur  and  Pollard 
(1986)  employed  a  more  realistic  assumption,  that  of  the  solution-precipitation  path. 
To  sustain  the  charge  transfer,  reactants  must  be  brought  to  the  reaction  layer  and 
the  products  removed.  As  the  charge  transferred  locally  increases  and  the  solubility 
limits  are  exceeded,  the  nucleation  process  commences  and  proceeds  with  the  rate 
given  by 

r„  =  K[{crC2)  -  IQ  (12) 

where  kn  is  the  rate  constant  for  the  nucleation  process,  (ciCi)  represents  the  local 
supersaturation,  and  Kg  is  the  solubility  product.  With  further  passage  of  current, 
the  growth  of  crystallites  and  formation  of  new  nuclei  occur  simultaneously.  The  rate 
of  crystallite  growth  obeys 


rg  =  kgvQiciCi)  -  Kg]  (13) 

so  that  the  precipitation  rate  constant,  k,  for  the  formation  of  a  solid  phase  and,  ipso 
facto,  for  the  removal  of  Cl~  from  the  reaction  space,  is  given  by 

k  =  kn-kkgv"^.  (14) 


Solving  eq.  (14)  for  v  and  substituting  into  eq.  (11),  yields 


s{l,t)  =  ^ 


1  - 


S(/.0 

(15) 


The  term  u”*,  in  eqs.  (13)  and  (14),  describes  in  general  terms  the  crystal  sur¬ 
face  available  for  growth,  as  does  the  parameter  p{l,t)  in  eq.  (11),  the  exponent  m 
also  reflects  the  morphology  of  the  growing  crystallites:  large  values  of  m  character¬ 
ize  crystallites  with  a  small  surface-to-volume  ratio  and  vice  versa.  Thus,  spherical 
particles  would  have  a  large  m  while  needle-  and  plate-like  growths  would  be  charac¬ 
terized  by  small  m.  Equation  (15)  indicates  that  the  rate  of  nucleation  and  crystallite 
growth,  together  with  the  morphology  of  the  precipitate,  strongly  affect  the  cell  life¬ 
time.  Calculations  carried  out  for  the  LifSOCl^  cells  imply  that  the  cell  lifetime  can 
be  extended  by  fast  precipitation  of  LiCl.  Conceivably,  in  the  limit,  the  Li/SOCl2 
cell  lifetime  could  be  doubled  (Tsaur  k  Pollard,  1986). 
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The  importance  of  the  m  exponent  cannot  be  overstated.  For  a  specific  positive 
electrode  structure,  an  increase  in  jo  tends  to  decrease  cell  lifetime.  However,  an 
opposite  effect  was  found  upon  the  addition  of  transition  metal  phthalocyanines  and 
other  additives.  While  the  catalyst  certainly  affects  the  exchange  current  density, 
its  role  is  much  more  complex  since  cell  modeling  and  experimental  evidence  suggest 
that  the  extension  of  the  cell  lifetime  is  connected  with  a  change  in  the  morphology 
of  the  emerging  new  solid  phase  (Nowak  et  al.,  1988). 

3.3.1.  Model  Predictions  vs.  Experiment 

Experimental  and  computed  galvanostatic  discharge  curves  of  a  Li/SOCl2  cell 
employing  an  acid  electrolyte  are  shown  in  figure  20.  The  hump,  at  the  beginning  of 
the  discharge,  is  attributed  to  SOCI2  reduction  in  an  acid  electrolyte,  i.e.,  the  width 
of  this  hump  is  proportional  to  the  amount  of  excess  AICI3  initially  present  in  the 
discharging  cell  (for  a  detailed  discussion  of  the  mechanism  refer  to  section  2.6).  The 
slope  of  the  central  region  of  the  Uc{i,t)  curve  depends  primarily  on  the  properties 
of  the  LiCl  film  at  the  negative  electrode  surface.  Toward  the  end  of  the  discharge, 
after  approximately  65  percent  of  the  available  material  has  been  used,  the  shape 
of  the  Uc{i,t)  curve  is  governed  by  the  events  occurring  within  the  porous  structure 
of  the  positive  electrode.  In  agreement  with  observation,  neither  the  slope  nor  the 
cell  voltage  are  sensitive  to  the  sjo  product;  however,  this  does  not  apply  to  the  cell 
lifetime,  viz.,  the  larger  the  sjo  product,  the  shorter  the  lifetime.  The  effect  of  the 
sjo  product  on  the  cell  lifetime  can  be  examined  with  the  aid  of  an  electric  circuit 
analog  (Szpak,  1991). 


3. 3. 1.1.  Transfer  Current  Density.  In  the  presence  of  excess  AICI3,  the  calculated 
charge  transfer  current  density  profile  is  dictated  by  a  balance  between  the  effect  of 
the  electrode  matrix  to  solution  conductivities  ratio  and  the  sjo  product.  For  the 
reduction  of  SOCI2,  with  jo  =  fix  10~^  A  cm~^  (Madou  &  Szpak,  1984)  the  sjo 
product  is  low  enough  to  force  the  reaction  zone  deeper  into  the  electrode  porous 
structure  (Tsaur  &  Pollard,  1986;  Szpak,  1991).  The  evolution  of  the  current  density 
(cd)  distribution  with  time  is  illustrated  in  figure  21.  Initially,  the  distribution  is 
normal  and  characteristic  of  low  matrix  resistance  (Szpak  &  Driscoll,  1983).  With 
the  passage  of  time,  a  change  in  the  reaction  path  occurs  due  to  the  depletion  of 
AICI3  near  the  front  face  while  in  the  back  of  the  electrode  the  solution  remains 
acidic.  The  result  is  the  appearance  of  a  minimum  in  the  cd  curves,  curves  2.  3, 


Figure  20.  Schematic  representation  of  a  galvanostatic  discharge  curve;  effect  of  change 
in  kinetic  parameters  (Driscoll  et  al.,  1985). 


etc.,  which  becomes  more  pronounced  and  penetrates  deeper  into  the  structure  with 
further  passage  of  time. 


Reaction  rate  distribution  in  porous  electrodes  cannot  be  measured  directly.  Thus, 
one  can  only  verify  the  calculated  profiles  indirectly,  e.g.,  by  examining  the  dis¬ 
tribution  of  precipitated  LiCl  within  the  cathode  structure.  Electron  dispersive 
x-ray  analysis  (  EDAX  )  of  electrodes  from  partially  discharged  cells  employing  acid 
electrolyte  was  performed  elsewhere  (Szpak  h  Driscoll,  1983;  Driscoll  et  al.,  1985). 
These  results  did  show  that  the  charge  transfer  reaction  occurs  initially  in  the  elec¬ 
trode  segment  close  to  the  electrode/separator  interface,  as  suggested  by  the  model. 
Experimental  profiles  at  later  times,  show  the  development  of  a  weak  minimum  which, 
however,  does  not  penetrate  as  deep  as  calculated.  The  likely  reason  for  this  discrep¬ 
ancy  is  the  assumed  relation  between  active  surface  area  and  the  porosity. 


3. 3. 1.2.  Optimum  Electrolyte  Composition.  Transport  properties  in  the  elec¬ 
trolyte  phase  are  the  most  important  input  data  in  calculating  the  electrolyte  concen¬ 
tration  needed  to  optimize  the  cell  lifetime.  Excellent  agreement  between  the  model 
and  experimental  data  is  obtained  for  the  relationship  between  the  initial  electrolyte 
concentration  and  the  cell  lifetime.  As  shown  in  figure  22,  the  optimum  cell  lifetime 
for  neutral  solutions  is  predicted  to  be  at  a  concentration  of  1.5  M  LiAlCU  while  lab- 
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Figure  21.  Reaction  distribution  in  the  positive  electrode  at  different  utilizations: 
(1)  0.6%,  (2)  1.49%,  (3)  1.79%,  (6)  11.36%  (Tsaur  &  Pollard,  19846). 


oratory  data  show  that  1.6  M  is  the  optimum  concentration.  This  excellent  agreement 
has  been  realized  only  after  the  relevant  data  were  obtained  (Szpak  &  Venkatasetty, 
1983). 

3. 3. 1.3.  Cathode  Thickness.  On  the  basis  of  a  large  number  of  experiments, 
Driscoll  (1984)  established  that  the  optimum  power  densities  are  obtained  with  cath¬ 
ode  thickness  between  2.6  x  10““*  and  4.0  x  10““*  m  .  The  calculated  effect  of  the 
electrode  thickness  on  the  cell  lifetime  for  discharge  rates  comparable  to  the  experi¬ 
mental  conditions,  shown  in  figure  23,  is  in  excellent  agreement.  Evidently,  increasing 
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Figure  22.  Calculated  dependence  of  cell  lifetime  on  electrolyte  composition  (Tsaur  & 
Pollard,  1984a). 


the  thickness  beyond  the  value  of  Lp2  =  4.0  x  lO”"*  m  does  not  improve  cell  life¬ 
time  (or  power  output)  significantly  because  the  reaction  does  not  penetrate  the 
added  electrode  volume.  Such  increases  would  only  reduce  the  volumetric  energy 
and  power  density  of  the  cell.  On  the  other  hand,  reducing  the  thickness  below 
Lpi  =  2.6  X  10“^  m  lowers  the  cell  lifetime  dramatically  and  imposes  stringent 

quality  control  requirements  because,  in  this  region,  the  lifetime  is  sensitive  to  small 
changes  in  electrode  thickness. 


3. 3. 1.4.  Electrode  Swelling.  Expansion  of  the  positive  electrode  during  cell  dis¬ 
charge  was  first  reported  by  Marincic  in  a  series  of  papers  (Marincic,  Epstein,  k 
Goebel,  1975;  Marincic,  1976a,  19766;  Marincic  k  Koffmann,  1976).  Somewhat  later, 
Szpak  and  Driscoll  (1983)  re-examined  this  problem.  They  emphasized  the  impor- 
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Figure  23.  Calculated  dependence  of  cell  lifetime  as  a  function  of  electrode  thickness 
(Driscoll  et  al.,  1985). 


tance  of  porosity  and,  in  particular,  pore  size  distribution  and  concluded  that  a 
30-percent  increase  in  cell  lifetime  is  possible  if  space  for  electrode  expansion  during 
discharge  is  provided.  The  calculated  increase  in  cell  lifetime  is  much  higher,  viz.,  75 
percent,  provided  that  all  free  space  generated  in  the  course  of  cell  discharge  can  be 
utilized. 


3.4.  THERMAL  MANAGEMENT 


In  its  simplest  form,  thermal  management  requires  maintaining  the  battery  tem¬ 
perature  within  prescribed  limits  during  the  course  of  discharge.  Although  these 
limits  may  vary,  the  upper  limit  for  a  Li  battery  should  not  exceed  the  melting  point 
of  lithium  if  the  initiation  of  highly  exothermic  reactions  is  to  be  avoided.  However, 
the  temperature  range  has  not  been  firmly  established  as  yet.  Reportedly,  experience 
with  high-rate  batteries  indicates  that  410  K  is  the  highest  safe-operating  tempera¬ 
ture,  however,  claims  have  been  made  that  the  Li/SOCl2  interface  is  stable  above 
the  melting  point  of  lithium  (Babai,  Bababigi,  &:  Bineth,  1984). 
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The  rate  of  heat  generation  in  a  discharging  battery,  derived  from  the  first  and 
second  laws  of  thermodynamics,  is 

q=j[U{Q)-U{j)]  +  a  (16) 

where  cr  is  a  function  representing  the  rate  of  heat  generation  arising  from  other  than 
purely  kinetic  processes.  Thus, 


where  the  first  term  is  qent  for  the  m  charge  transfer  processes  and  the  second  term 
is  the  change  in  enthalpy  arising  from  the  n  chemical  reactions.  This  function  can 
ba  evaluated  experimentally  and  the  nature  and  the  magnitude  of  some  of  the  heat 
sources  assessed  (Godshall  &  Driscoll,  1984).  For  example,  it  has  been  shown  that 
localized  heat  source  due  to  corrosion  of  Li  in  SOCI2  increases  with  an  increase  in 
the  discharge  current  (Szpak  &  Driscoll,  1987). 


Since  the  total  heat  generation  of  a  discharging  battery  is  a  function  of  local 
current  density  and  temperature,  eqs.  (16)  and  (17),  it  follows  that,  to  obtain  the 
temperature  profiles  within  the  interior  of  an  operating  battery,  eqs.  (18)  and  (19) 
must  be  solved  simultaneously 


XVU  +  j  =  0  (18) 

dT 

=  div(KVr)  +  a.  (19) 

Equations  (18)  and  (19)  are  coupled  by  the  dependence  of  A  on  temperature  and 
of  cr  on  both  temperature  and  solution  conductivity.  In  practice,  it  is  convenient  to 
express  the  heat  generation  in  terms  of  discrete  sources,  i.e., 

n 

with  the  main  source  being  attributed  to  the  joule  heating  in  the  cell  functional 
elements,  i.e.,  polarization  of  individual  electrodes  and  internal  cell  resistance.  On 
occasion,  additional  heat  sources  may  be  identified,  e.g.,  due  to  intercell  currents 
always  present  in  the  reserve  type  design  and  chemical  reactions  between  the  cell 
components  and  reaction  products. 
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The  boundary  and  initial  conditions  must  reflect  the  mode  of  battery  discharge 
and  be  consistent  with  battery  design.  Parnell  and  Szpak  (1985)  reported  on  the 
temperature  distribution  within  a  module  consisting  of  a  large  number  of  identical, 
disc-shaped  cells  of  bipolar  construction.  In  such  a  design,  heat  transport  in  the  axial 
direction  is  minimal,  except  at  the  terminal  cells.  In  practice,  this  leaves  only  one 
spatial  dimension  to  consider. 

3.4.1.  Model  Predictions  vs.  Experiment 


In  contrast  to  batteries  containing  aqueous  electrolytes  where  solutions  of  eqs.  (18) 
and  (19)  with  constant  coefficients  yield  satisfactory  results,  realistic  analysis  of  the 
Li/SOCli  battery  requires  that  these  coefficients  reflect  changes  that  have  occurred 
in  the  electrolyte  phase,  in  particular,  changes  in  solution  composition  and  tempera¬ 
ture.  The  use  of  average  values  for  thermal  and  electrical  properties  results  in  higher 
cell  temperatures  and  larger  temperature  gradients.  The  practical  implications  are 
obvious:  uniform  temperature  profile  indicates  uniform  cd  distribution  and,  conse¬ 
quently,  uniform  consumption  of  reactants  and  uniform  reduction  in  the  porosity  of 
the  positive  electrode.  In  addition,  because  of  lower  cell  temperature,  either  higher 
cell  currents  or,  alternatively,  longer  discharge  times  can  be  realized  before  the  cell 
interior  temperature  exceeds  an  a  priori  set  limit.  Parnell  and  Szpak  (1985)  suggested 
that  to  obtain  agreement  with  experiment,  the  model  must  consider  changes  in  the 
individual  elementary  processes  that  take  place  during  the  cell  discharge,  i.e.,  changes 
in  the  coefficients  and  sources  in  eqs.  (18)  and  (19).  This  inclusion  renders  the  differ¬ 
ential  equations  nonlinear  and  forces  a  numerical  solution  of  the  relevant  equations. 
Using  this  procedure,  good  agreement  between  the  computed  and  observed  tempera¬ 
ture  hcis  been  obtained,  figure  24.  They  also  found  that  in  the  cathode  limited  cells 
the  prime  heat  source  resides  in  the  cathode  structure  when  the  cathode  porosity 
becomes  less  than  ca  20  percent. 


3. 4. 1.1.  Activation  Under  Load.  As  a  rule,  a  high-rate,  reserve  battery  is  acti¬ 
vated  at  open  circuit.  There  are,  however,  operational  situations  requiring  activation 
under  load.  To  simulate  this  requirement,  the  rate  of  electrolyte  penetration  into  the 
cell  must  be  specified.  Parnell  and  Szpak  (1985)  expressed  the  penetration  rate  in 
the  form  of  the  time-dependent  increase  in  the  electrode  active  area  by  an  arbitrarily 
selected  expression, 

«(0  =  «max(l-e-''‘)  (20) 
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Figure  24.  Computed  and  measured  temperatures  at  the  inner  cell  diameter: 
(a)  experimental,  (b)  computed  for  variable  coefficients  (Parnell  &  Szpak,  1985). 


where  the  0  parameter  determines  the  rate  at  which  the  active  segment  of  a  cell  grows 
or,  alternatively,  the  rate  at  which  local  cd  decreases  with  time. 

A  plot  of  temperature  versus  time  at  the  point  of  the  electrolyte  entry  is  shown  in 
figure  25.  Evidently,  the  slower  the  fill  rate,  the  higher  the  initial  temperature.  When 
the  cell  is  completely  filled,  the  effect  is  less  and  is  indistinguishable  after  about  6 
minutes.  The  major  conclusion  concerning  safety  is  that  activation  under  load  should 
be  discouraged  unless  precautions  are  taken  to  cissure  a  fast  and  reliable  fill. 


CELL/BATTERY  MODELING  -  CATASTROPHIC  THERMAL 
RUNAWAY 


In  this  system,  designed  for  high  rates  of  discharge,  violent  reactions  have  occurred 
during  the  development  process.  A  number  of  explanations  have  been  offered,  but 
with  very  little  substantiating  data  except  for  clear  cases  of  thermal  runaway.  Szpak, 
Gabriel,  and  Driscoll  (1987)  examined  the  initiation  and  progress  of  the  catastrophic 
thermal  excursion  by  considering  a  heat  source  extrinsic  to  the  normal  mode  of  cell 
operation,  i.e.,  heat  source  q\{r)  representing  joule  heating  arising  from  an  internal 
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Figure  25.  Temperature/time  relationship  as  a  function  of  filLng  rate  under  load 
(Parnell  &  Szpak,  1985). 


short  circuit.  As  the  temperature  rises  locally  above  the  critical  temperature  T^, 
this  source  activates  the  intrinsic  ,  associated  with  the  exothermic  reaction 

involving  lithium  and  other  cell  components  e.g.  ,  SOCh 

92  =  ~-Mlrkocexp{-E/RT)^T  -  T,,).  (21) 


Equation  (21)  expresses  the  strength  of  the  intrinsic  heat  source  as  a  product  of 
the  heat  of  reaction,  —AHr-,  and  the  reaction  rate.  The  Heaviside  unit  step  function 
is  employed  to  account  for  the  protective  action  of  the  LiCl  film,  always  present 
whenever  metallic  Li  is  in  contact  with  SOCI2  ,  i.e., 

$(T  -  7;,)  =  0;  for  T  <  Ter 


$(T  -  Ter)  =  1;  for  T  >  T^- 
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The  concentration  term  in  eq.  (21)  must  satisfy  the  diffusion  equation 


—  =  div(£)Vc)  -  k^cexY>{-EIRT) 


(22) 


where  the  sink  term  is  the  reaction  rate.  To  complete  the  mathematical  description, 
the  initial  and  boundary  conditions  must  be  specified.  They  assume  as  initial  con¬ 
ditions  uniform  temperature  and  solution  concentration  and  the  boundary  condition 
at  r  =  a  are  the  fluxes  given  by  eqs.  (23)  and  (24) 

dT 

-\—  =  h,{T-To),  (23) 


(24) 


These  conditions  allow  heat  transport  in  the  radial  direction  but  disallow  mass  trans¬ 
port  across  this  boundary. 


Evolution  of  temperature  profiles  with  time  following  the  activation  of  the  extrin¬ 
sic  heat  source  are  shown  in  figure  26.  At  first,  the  temperature  profiles  are  those 
associated  with  centrally  localized  heat  sources,  curves  1  and  2.  As  soon  as  the  criti¬ 
cal  temperature  is  exceeded,  the  temperature  profile  changes  to  that  arising  from  two 
sources:  9i(r)  and  the  time-dependent  intrinsic  source  q2{r',t)  causing  the  distortion 
in  temperature  profiles,  curves  3  and  4.  With  a  further  passage  of  time,  the  tempera¬ 
ture  profiles  take  on  a  skewed,  spike-like  appearance,  indicating  the  narrowing  of  the 
reaction  zone  and  a  rapid  increase  in  the  heat  intensity,  curves  5  and  6,  followed  by 
a  period  of  stabilization,  curves  7  and  8. 


The  catastrophic  event  is  characterized  by  a  period  of  induction,  followed  by  ig¬ 
nition  and  propagation  of  a  narrow,  very  hot  zone  of  burning  Li  with  temperatures 
well  above  the  melting  points  of  aluminum  and  stainless  steel.  An  example  of  ther¬ 
mal  runaway  caused  intentionally  in  a  2000-.Mi  cell  is  shown  in  six  photographs  in 
figure  27.  An  analysis  of  numereous  videotapes  of  battery  ventings  was  consistent  with 
these  features:  the  existence  of  an  induction  period,  followed  by  ignition,  propagation 
of  a  reaction  zone  and  flame  extinction  which,  as  a  rule,  follows  the  venting  of  a  cell. 
Furthermore,  the  evolution  of  temperature  profiles  using  liquid  crystals  deposited  on 
an  external  cell  surface  follows  such  a  sequence  (Driscoll  Szpak,  1985). 
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Figure  26.  Evolution  of  temperature  (solid  lines)  and  concentration  (dashed  lines) 
following  the  activation  of  an  extrinsic  heat  source  (Szpak  et  al.,  1987). 


3.6.  MODULE/BATTERY  MODELING  -  INTERCELL  CURRENTS 


By  definition,  an  intercel!  current  is  an  ionic  current  that  originates  in  one  cell 
and  terminates  in  another  (Szpak,  Gabriel,  &  Driscrll,  1984).  This  current  does  not 
perform  any  useful  work;  on  the  contrary,  it  is  a  parasitic  current  resulting  not  only 
in  a  loss  of  power  but  also  in  a  loss  of  capacity,  and  may  contribute  to  the  premature 
failure  of  individual  cells  by  the  formation  of  passive  films,  dendritic  growths,  or  by 
the  corrosion  of  structural  parts. 


In  a  discharging  cell,  two  processes  occur  at  the  negative  electrode:  charge  trans¬ 
fer  and  transport  of  Li'^  ions  through  the  protective  LiCl  film.  The  thickness  of 
this  film  depends  on  the  electrolyte  composition  and  the  rate  of  battery  discharge. 
While  the  situation  at  the  positive  electrode  is  far  more  complex,  we  assume  that 
each  elementary  process  obeys  a  generalized  Ohm’s  law,  i.e.,  it  can  be  described  by  an 
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equation  of  the  type,  j  =  Ary,  relating  the  flux  j  to  the  driving  force,  y,  through  a 
phenomenological  rate  constant,  k.  Such  a  description  permits  the  construction  of  an 
electric  circuit  analog,  shown  in  figure  28,  to  represent  the  overall  process  occurring 
in  the  cell  interior.  It  is  noteworthy  that  this  type  of  representation  has  served  as  the 
starting  point  in  practically  all  cases  reported  in  the  literature  (Ksenzhek  &  Koshel, 
1970,  1971;  Nesterov,  Razevig,  &  Korovin,  1974;  Onishchuk,  1972;  Thiele,  ScheifF, 
&  Matschiner,  1983;  Kaminski  &;  Savinell,  1983).  This  approach  provides  a  simple 
means  of  treating  intercell  currents  that  yields  results  in  agreement  with  observation. 


R| 


Figure  28.  Battery  representation  by  equivalent  electric  circuit  analog.  cell 

voltage;  potential  at  manifold  inlet  to  j  —  1-th  cell;  J,:  circulating  current; 

J/:  loop  cnrrent  in  an  external  resistor  /?/;  Rj:  manifold  feedline  resistor;  Rf  .  equiv¬ 
alent  feed-tube  resistor.  Control  volume  indicated  by  dashed  line  (Szpak  et  al.,  1984). 


Assuming  the  constancy  of  phenomenological  coefficients  associated  with  the  ionic 
current  and  replacing  them  with  equivalent  resistances,  /?<,,  and  /?/,,.  the  loop  equa¬ 
tions  for  the  current,  .7;  and  an  external  load,  /?;,  are  given  by 

=  =  l,2,.3,...,yV.  (2.')) 
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The  circulating  currents  in  the  fill  path,  J,,  obey 
—  +  —  —  =  1, 1, iV— 1.  (26) 


Additional  statements  needed  to  complete  the  circuit  description  are:  circulating 
currents  in  the  zero-th  and  N-th  segments,  Jo  =  Jn  =  0;  the  bipolar  plate  is  equipo- 
tential,  i.e.,  K.i+i  =  K;,t;  the  total  anodic  current,  e.g.,  in  the  i-th  cell  is  a  sum  of 
the  load  current  Ji  and  the  circulating  current  J<+i;  cind  the  intercell  current  /,  is  the 
difference  in  the  circulating  currents  in  the  i-th  and  i-1  th  segments,  li  =  Ji  —  Ji-i- 


Equations  (25)  and  (26)  were  formulated  for  a  discrete  model.  They  can  be 
reformulated  for  a  continuous  model  that  has  the  advantage  of  admitting  solutions 
in  a  closed  form  (Ksenzhek  &  Koshel,  1970,  1971).  The  selection  between  these 
approaches  depends  upon  the  desired  accuracy  and  the  specifics  of  the  electrochemical 
system.  For  example,  the  continuous  model  may  not  be  relevant  for  calculating 
intercell  currents  in  batteries  where  operation  terminates  due  to  blocked  electrode 
surface  (Szpak  et  al.,  1984).  For  the  simple  case,  that  with  constant  coefficients, 
solutions  to  eqs.  (25)  and  (26)  can  be  obtained  in  a  closed  form,  viz., 

N{l-i)U 
‘  Ri  +  Nil  -  ORz 

(  -f  A*  \ 

<7i  =  (  1  J  ]  )  *  =  1>  2, ...,  N 


(27) 

(28) 


where 


(29) 


A  =  1  +  r  +  s/n2  +  V)  ;  y  = 


K,  r  _  (A'"-i)(A  +  i)  • 
«,  +  «=[  A'(A'''  +  1)(A-1) 
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and 


fl;  =  R,  +  fl,  -  H.  (^) 


For  a  given  number  of  cells,  the  relative  distribution  of  both  the  circulatii  •  'd 
intercell  currents  depends  on  the  value  of  the  parameter  A.  This  parameter  is  a 
function  of  the  ratio  (/?/  +  R^)/ i.e.,  it  depends  on  the  battery  design  as  well  as 
on  the  properties  of  the  electrolyte  and  the  kinetics  of  the  charge  transfer  reaction. 
In  particular,  for  large  values  of  /?(,  A  will  be  large  while  for  large  values  of  /?/,  A  will 
be  near  unity.  In  the  latter  case,  the  tails  of  the  intercell  current  distributions  extend 
into  the  battery  with  a  shape  depending  on  the  number  of  cells.  In  contrast,  for  large 
A  ,  the  terminal  cells  carry  most  of  the  intercell  currents  independently  of  the  number 
of  cells  in  a  module.  In  practice,  it  is  desirable  to  have  large  Rt  because  it  minimizes 
the  circulating  currents,  J,,  which  means  the  intercell  currents  in  the  terminal  cells 
are  also  small  and  the  current  reversal  does  not  occur,  i.e.,  the  dendritic  growth  and 
the  corrosion  of  bipolar  plates  is  avoided. 


Often,  to  meet  performance  requirements,  a  multimodule  battery  is  assembled. 
An  example  of  a  multimodulc  assembly  is  shown  in  figure  29.  In  general,  an  assembly 
can  contain  either  an  odd  or  an  even  number  of  modules  which,  when  connected  in 
parallel,  share  positive/positive,  (+/+),  and  negative/negative,  (  —  /  —  ),  interfaces.  In 
an  assembly  of  an  odd  number  of  modules,  only  one  configuration  is  possible,  viz.,  one 
positive  and  one  negative  assembly  terminal.  For  an  even  number  of  modules,  two 
configurations  are  possible  both  assembly  terminals  are  either  positive  or  negative. 


The  formulation  of  governing  equations  is  the  same  as  for  the  single  module  except 
for  the  matching  conditions  at  the  (+/+)  and  the  (  — /  — )  interfaces  between  the 
neighboring  modules.  To  avoid  ambiguity,  quantities  associated  with  a  module  are 
identified  here  by  a  preceding  index  in  parenthesis,  e.g.,  *"’1 J/  denotes  the  load  current 
generated  by  the  m—  th  module.  For  an  A/—  module  battery,  the  potential  across 
the  load  is, 

M 

V  ^  RiJ2  (30) 

m=l 

As  before,  the  circvilating  currents  in  terminal  cells  of  a  battery  are  zero.  However,  in 
terminal  cells  of  a  module,  i.e.,  at  coupling  interfaces  between  modules,  the  circulating 
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Figure  29.  Schematic  representation  of  a  multimodule  battery,  (a)  solid  lines  for  an 
odd  number  of  modules,  (b)  dashed  lines  for  an  even  number  of  modules  (Szpak  et  al., 
1990). 


currents  are  not  necessarily  zero;  in  particular  they  are 

=  -I-Un  = 

"Jo  =  =  J, 

where  J(m-i),(m)  and  t/(Tn),(m+i)  are  the  coupling  currents  between  modules.  The 
mathematics  is  straightforward  but  somewhat  lengthy.  The  reader  is  referred  to 
original  work  (Nesterov  et.  al.,  1974). 


3.6.1.  Model  Predictions  vs.  Experimental  Observation 


A  number  of  specific  situations  can  be  e.xamined.  Here,  we  limit  the  discussion 
to  the  distribution  of  circulating  and  intercell  currents,  conditions  leading  to  current 
reversal,  and  reduction  of  power  output.  We  complete  the  discussion  by  examining 
the  behavior  of  a  multimodule  battery. 

3. 6. 1.1.  Distribution  of  Intercell  and  Circulating  Currents.  The  distribution 
of  intercell  and  circulating  currents  in  a  LifSOCl2  multicell  module  is  shown  in  figure 
30.  It  is  seen  that,  independent  of  the  number  of  cells  in  the  module,  maximum 
circulating  currents  are  at  the  central  cells,  whereas  the  intercell  currents  are  largest 
at  the  ends.  Also,  the  intercell  currents  change  their  signs  from  positive  to  negative, 
i.e.,  from  entering  the  cell  at  the  negative  end  of  a  module  to  exiting  the  cell  at 
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the  positive  ends.  Evidently,  in  a  module  having  a  small  number  of  cells,  all  cells 
generate  intercell  currents.  However,  as  the  number  of  cells  increases,  the  centrally 
located  cells  contribute  less  and,  in  the  case  of  an  80-cell  module,  the  majority  of  cells 
contribute  nothing  at  all  to  the  intercell  current.  The  practical  implication  is  that, 
as  the  number  of  cells  increases,  losses  in  the  extractable  energy  approach  a  limiting 
value  and  that  most  of  the  module  suffers  little,  if  any,  from  reactant  underutilization. 
The  shape  of  the  distribution  curve  of  the  intercell  current  implies  that  a  battery  may 
be  viewed  as  consisting  of  segments,  viz.,  a  central  segment  where  cells  operate  with 
equal  anodic  and  cathodic  currents  and  segments  involving  terminal  cells  where  the 
underutilization  of  electroactive  materials  occurs. 


Figure  30.  Dist  ribution  of  circulating  (thin  line)  and  intercell  currents  (heavy  line)  as 
a  function  of  the  number  of  cells;  (a)  small  number  of  cells  {N  =  10),  (b)  moderate 
number  of  cells  (N  =  40),  (c)  large  number  of  cells  {N  =  80).  Current  reversal 
indicated  (Szpak  et  al.,  1984). 


3. 6. 1.2.  Cell  Current  Reversal.  Szpak  et  al.  (1984)  have  shown  that,  if  the  inter¬ 
cell  current  exceeds  a  critical  value,  Icr,  which,  for  a  given  set  of  design  parameters, 
depends  on  the  module  load,  the  cell  current  can  be  reversed  over  a  part  of  an  elec¬ 
trode  in  one  or  more  cells,  as  shown  in  figure  30c.  In  particular,  for  the  current  into 
the  cell,  the  cell  current  reversal  occurs  on  the  anode,  while  for  the  exiting  intercell 
current,  it  takes  place  on  the  cathode.  Thus,  at  the  negative  end  of  a  module,  in  a 


number  of  cells,  dendritic  growth  of  metallic  Li  is  possible  and  is  frequently  observed, 
while  the  positive  electrode,  in  a  corresponding,  but  not  necessarily  equal  number 
of  cells,  suffers  either  corrosion  or  passivation.  The  effect  of  cell  current  reversal  is 
shown  in  figure  31  where  photographs  taken  in  the  vicinity  of  the  feedline  region  show 
distinct  evidence  of  electrochemical  activities  at  the  positive  and  negative  ends  of  a 
module  with  no  measurable  changes  in  the  central  part. 


Figure  31.  Detrimental  effect  of  intercell  current  at  the  cell  entrance  and  within  the 
annular  feedline;  (top  left)  cell  located  at  the  negative  end  -  dendritic  growth  of 
metallic  lithium  indicated  by  arrow,  (top  tight)  cell  located  in  the  central  region  of 
an  80-cell  module  -  no  visible  changes,  (bottom)  cell  located  at  the  positive  end  - 
corrosive  attack  indicated  by  an  arrow  (Szpak  et  al.,  1984). 


3. 6. 1.3.  Power  Losses.  The  presence  of  a  common  electrolytic  path  also  adversely 
affects  the  module  power  output.  Gabriel  and  Szpak  (1989)  examined  power  losses 
in  a  simple  case  of  a  module  containing  N  identical  cells  represented  by  an  electric 
circuit  analog  with  constant  coefficients.  .A,  plot  of  dimensionless  power,  Ilo,  as  a 
function  of  dimensionless  resistance,  p,  shown  in  figure  32,  indicates  that  the  power 
into  the  load  exhibits  a  maximum  at  p  =  1  —  C,  while  battery  own  power  (Godshall 
&  Driscoll,  1984)  decreases  monotonically  with  increasing  p.  It  is  noted  that  for  a 
given  the  power  consumption,  fib  —  flj  increases  from  (  to  unity  with  decreasing 
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p.  This  simple  model  accounts  for  dissipative  losses  occurring  during  the  discharge 
and  gives  their  distribution  for  small  intercell  currents.  The  initial  distribution  of 
dissipative  processes  is  symmetrical  with  a  shift  in  the  distribution  occurring  in  the 
course  of  module  discharge.  The  magnitude  of  this  shift  depends  on  the  charge 
transfer  reactions  and  the  manifold  design.  Defective  cells  may  create  highly  localized 
heat  sources  which,  in  turn,  may  initiate  a  catastrophic  event. 


Figure  32.  Plot  of  dimensionless  power,  Ilo,  with  a  =  f,  6  as  a  function  of  dimensionless 
resistance,  p,  for  ^  =  0.0,  0.1,  and  0.2,  as  indicated  (Gabriel  &  Szpak,  1989). 


A  technologically  important  aspect  of  battery  scale-up  is  its  effect  on  the  distribu¬ 
tion  of  load  currents  from  individual  modules.  The  addition  of  modules  reduces  the 
individual  load  currents  while  maintaining  a  smooth  current  delivery  to  the  main  line, 
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except  toward  the  end  of  battery  discharge  when  crossovers  among  module  currents 
occur,  indicating  current  sharing  that  does  not  interfere  with  the  power  output.  As 
illustrated  in  figure  33,  the  measured  delivery  of  total  current  is  smooth,  but  the  cur¬ 
rent  supplied  by  individual  modules  varies.  The  experimentally  measured  currents  in 
a  three-module  battery,  shown  in  figure  34,  exhibit  features  that  are  in  agreement  with 
calculated  individual  and  nodal  currents.  The  addition  of  modules  to  a  well-designed 
and  constructed  multimodule  battery  does  not  appear  to  increase  safety  concerns, 
i.e.,  does  not  inherently  cause  an  increcise  in  the  imbalance  of  intercell  currents.  By 
well-designed  and  constructed  batteries,  we  mean  the  following  restrictive  conditions 
apply:  all  modules  contain  the  same  number  of  cells  and  the  construction  of  the  cells 
and  their  performances  are  statistically  reproducible. 


Figure  33.  Calculated  load  current  variation  for  three-module  battery  assembly  (for 
input  parameters  refer  to  Gabriel,  et  al.,  1990). 
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(a) 


(b) 


MINUTES 

Figure  34.  Nodal  load  current  variation  in  an  experimental  three-module  assembly; 
(a)  electric  circuit  analog  of  a  three-module  battery,  (b)  measured  current  variation 
during  battery  discharge  (Szpak  et  al.,  1990). 
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4.  FABRICATION  AND  QUALITY  CONTROL 


4.1.  mTRODUCTORY  COMMENTS 


The  problems  encountered  in  fabrication  and  quality  control  become  more  diffi¬ 
cult  to  resolve  as  the  expected  performance  level  is  increased.  To  meet  the  design 
criteria,  i.e.,  to  maximize  the  extractable  energy  content  for  a  given  power  required,  a 
set  of  optimizing  procedures  is  performed.  These  optimization  procedures  favored  a 
thin  cell  design.  In  practice,  this  means  separating  all  cell  components  with  minimum 
ohmic  losses.  The  efficient  operation  of  reactors  designed  for  electrochemical  process¬ 
ing,  including  energy-conversion  devices,  often  requires  the  installation  of  a  common 
electrolytic  path.  There  are  numerous  advantages  associated  with  this  approach,  viz., 
bipolar  cell  construction,  controlled  supply  of  reactants  and  removal  of  products,  as 
well  as,  a  reliable  means  of  thermal  management.  These  advantages,  however,  are 
partially  offset  by  the  parasitic  shunting  action  of  the  intercell  currents  (Szpak  et  ah, 
1984).  Briefly,  these  currents  not  only  reduce  the  power  output  but  often  affect  the 
structural  integrity  of  battery  functional  elements.  Modeling  and  practical  experi¬ 
ence  have  shown,  that  in  the  case  of  well-designed  modules  connected  electrically  in 
parallel  and  hydraulically  in  series,  an  increase  in  the  number  of  modules  does  not 
create  serious  problems  (Szpak  et  ah,  1990).  However,  in  an  assembly  of  modules 
with  one  or  more  defective  cells,  or  with  a  sufficiently  large  variation  among  the  cell 
functional  elements,  a  situation  arises  that  substantially  reduces  the  battery  lifetime 
and,  on  occasion,  may  lead  to  catastrophic  events. 


4.2.  MODULE  ASSEMBLY/QC  POINTS 


From  the  point  of  view  of  battery  fabrication  and  operation,  two  interrelated 
areas  require  special  attention.  The  first  area  involves  problems  connected  with  the 
activation  time  for  which  an  internal  distribution  system  for  the  flow  of  electrolyte 
must  be  devised.  Usually,  more  than  one  fill  tube  running  through  the  length  of  the 
module  are  needed.  Short  activation  times  of  a  large  battery  present  special  difficul¬ 
ties  because,  to  realize  uniform  and  fast  loading,  two  sets  of  mutually  contradicting 
requirements  must  be  simultaneously  satisfied.  First,  the  separator  must  be  porous 
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enough  to  permit  a  rapid  flow  of  electrolyte.  Yet,  the  pore  size  must  be  small  enough 
to  prevent  the  penetration  of  carbon  particles  and  eliminate  the  development  of  inter¬ 
nal  shorts.  In  addition,  the  separator  must  be  mechanically  strong  to  prevent  tearing 
and  shape  distortion  by  the  forced  flow  of  electrolyte  during  activation.  To  illustrate 
the  inherent  difficulties  caused  by  a  rapid  flow  of  electrolyte,  note  that  a  tenfold 
decrease  in  the  separator  porosity  necessitates  a  hundredfold  increase  in  pressure  to 
maintain  the  same  flow  rate. 


The  second  area  requiring  attention  is  thermal  management.  Often,  the  heat 
capacity  of  the  battery  is  insufficient  to  maintain  the  operating  temperature  within 
safe  limits  during  extended  periods  of  high  discharge.  Since  a  considerable  amount 
of  heat  must  be  rejected  to  the  surroundings,  either  a  highly  conductive  path  must 
be  provided  between  the  cell  interior  and  the  module  external  heat  transfer  surface, 
while  concurrently  maintaining  electrical  isolation  of  the  energy  section  from  the  heat 
rejection  surface,  or  a  flowing  electrolyte  must  be  used. 


4.3.  FABRICATION  OF  CELL  FUNCTIONAL  COMPONENTS 


A  typical  fabrication  sequence  of  a  module  and  the  pertinent  quality  control  (QC) 
points,  are  summarized  in  a  block  diagram  in  figure  35.  The  cell  functional  elements, 
e.g.,  positive  and  negative  electrode  structures,  are  first  fabricated.  Then,  other  com¬ 
ponents  such  as  separators,  bipolar  plates,  current  collectors,  etc.,  are  cut  to  size, 
followed  by  the  process  of  stacking  the  cells  and  compressing  the  stack  to  the  re¬ 
quired  thickness.  Finally,  a  means  for  battery  activation  or  electrolyte  circulation  are 
provided  to  complete  the  process  of  module  assembly.  The  compliance  with  specifi¬ 
cation  demands  counting  the  number  of  cells  in  a  module,  measuring  the  cell/module 
thickness,  and  inspecting  the  quality  of  seals.  Quality  control  measures  are  closely 
connected  with  the  fabrication  steps  and  the  dynamics  of  the  cell/module  operation 
(Gabriel  et  al.,  1990).  In  the  course  of  component  fabrication,  both  the  amount  of 
stored  electroactive  material  and  the  geometrical  constraints  of  the  porous  structure 
are  controlled.  The  control  of  the  purity  and  amount  of  the  electroactive  material 
presents  no  special  problems.  On  the  other  hand,  the  control  of  the  geometrical  con¬ 
straints  of  a  porous  electrode  matrix  is  a  serious  and,  in  many  instances,  unresolved 
problem. 
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FABRICATION  DISCHARGE 


Figure  35.  A  block  diagram  showing  the  fabrication  sequence  for  module/battery 
assembly  and  QC  requirements:  I  -  thickness,  N  -  number  of  cells,  Rf  and  Rt  - 
feedline  and  fill-tube  resistances  respectively,  S*  -  electrode  separation,  €  -  porosity 
(Gabriel  et  al.,  1990). 

4.3.1.  Positive  Electrode 

While  fabricating  the  positive  electrode,  the  following  factors  need  to  be  consid¬ 
ered:  large  specific  surface  area,  uniform  loading,  and  sufficient  mechanical  strength. 
Fabrication  methods,  and  the  resulting  properties  of  the  positive  electrode,  play  a 
crucial  role  in  achieving  the  required  performance  level.  The  fabrication  of  these 
electrodes  is,  as  a  rule,  protected  by  patent  rights  and  guarded  as  proprietary  infor¬ 
mation.  The  most  common  matrix  for  the  high-discharge  rate  Li/SOCl2  battery  is 
Teflon-bonded  Shawinigan  black  ( ontaining  various  additives  such  as  dispersed  cop¬ 
per  (Giattino  &  Catanzarite,  1979),  platinum  (Klinedinst,  1981),  and  iron  and  cobalt 
phthalocyanines  (Venkatasetty,  1981).  The  role  of  binders  is  to  assure  the  mechanical 
strength  while  the  intended  role  of  additives  is  to  modify  the  current  distribution,  by 
changing  either  the  reaction  path  or  the  growth  habits  of  the  precipitated  reaction 
product(s),  or  both. 

Quality  control  is  limited  to  a  strict  adherence  to  the  prescribed  fabrication  pro¬ 
cedures  and  the  quality  of  y?'\  materials  used.  Uniform  porosity  and  intimate  contact 
between  the  electrode  matrix  and  the  bipolar  plate  is  necessary  to  prevent  film  for¬ 
mation  by  elemental  sulfur  (Nowak  et  al..  1988). 
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4.3.2. 


Bipolar  Plate 


The  bipolar  plate  serves  two  functions:  it  provides  electronic  conductivity  between 
adjacent  cells  and  it  forms  a  thermal  path  for  rejecting  heat  from  the  cell  interior. 
The  construction  material  must  be  chemically  inert  with  respect  to  the  electrolyte. 
Nickel  and  stainless  steel  foils  satisfy  these  requirements.  However,  for  operating  at 
high-discharge  rates,  Ni—  clad  Cu  foils  are  used  to  facilitate  heat  transfer.  Quality 
control  procedures  include  visual  inspection  and  checks  for  specified  dimensions. 


4.3.3.  Lithium  Foil 


Battery  grade  metallic  lithium  foil  of  various  thicknesses  is  available  and  shipped  in 
steel  containers  under  an  argon  atmosphere.  Once  the  container  is  open,  the  material 
must  be  handled  in  a  dry  room  with  less  than  two  percent  relative  humidity. 

4.3.4.  Thionyl  Chloride  Electrolyte 

Thionyl  chloride  is  a  common  material  available  in  either  technical  or  purified 
grades.  As  a  rule,  thionyl  chloride  is  further  purified  by  distillation  over  Li  chips, 
dipentene,  or  triphenyl  phosphite.  It  is  noteworthy  that  impurities  have  a  negligible 
effect  on  the  performance  of  high-discharge  rate,  reserve  batteries.  However,  for  active 
cells,  impurities  can  and  do  have  a  deleterious  effect,  primarily  during  storage  (Driscoll 
et  al.,  1976).  The  electrolyte  is  prepared  by  dissolving  LiCl  and  AI2CIG,  which  are 
purchased  as  anhydrous,  reagent  grade  salts  and  are  normally  used  without  further 
purification. 


4.3.5.  Separator 


Separator  materials  include  glass,  ceramic,  and  organic  porous  sheets.  Tiie  func¬ 
tion  of  a  separator  in  a  battery  is  twofold:  to  assure  and  maintain  a  positive  separation 
between  the  negative  and  positive  electrodes  and  to  provide  an  ionically  conductive 
path.  In  a  specialty  batteiy  an  additional  requirement  is  imposed:  sufficient  porosity 
to  meet  activation  time  requirements  and  mechanical  strength  to  prevent  damage 
during  activation. 


4.4. 


MODULE  ASSEMBLY 


The  module  is  the  basic  building  block  in  the  design  of  a  multikilowatt  LilSOCl2 
battery  employing  a  static  electrolyte.  In  addition  to  cells  needed  to  produce  electric 
power  and  a  means  to  deliver  it  to  an  e.xternal  circuit,  the  module  also  contains  the 
hull  section.  The  basic  module  assembly  procedure  involves:  stacking  the  individual 
elements,  compressing  the  stack  to  the  final  dimension,  sealing  the  outer  and  inner 
cell  edges,  installing  the  hull  section  and  forming  the  thermal  joint.  The  fabrication 
sequence  is  summarized  in  figure  36.  An  endplate,  insulator  plate,  and  collector  plate 
are  placed  in  sequence  on  the  lower  platen  of  the  press.  This  is  followed  by  assembly  of 
the  cell  components,  viz.,  the  cathode,  separator,  and  an  anode/bipolar  subassembly. 
This  sequence  is  repeated  until  the  required  number  of  cells  is  assembled.  The  last 
anode  {Li  foil)  is  attached  to  the  collector  plate.  Installing  an  insulator  plate  and  an 
endplate  complete  the  module  stacking.  The  component  stack  is  then  compressed  to 
within  a  fraction  of  a  millimeter  of  the  final  height.  Bolts  and  spacers  are  inserted 
and  the  stack  is  compressed  to  the  final  height.  After  the  module  is  compressed  and 
bolted,  it  is  ready  for  sealing. 


INSTRUMENTED 
BI-POLAR  PLATE 

TEFLON  (MID-HEIGHT  OF  STACK) 


@  (EXTERNAL) 

Figure  36.  Exploded  view  of  a  module  section  illu.strating  the  installation  sequence  of 
cell  components.  Numerals  indicate  the  location  of  thermocouples. 
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DISTANCE 


Figure  38.  Method  preventing  internal  short  generated  by  intercell  current.  (A)  Li 
foil,  (B)  protective  LiCl  film,  (C)  porous  carbon  cathode,  (D)  bipolar  plate,  (E)  teflon 
disks,  (I)  intercell  current,  ia,  ic  ~  unodic  and  cathodic  cell  currents;  6^  -  effective  cell 
thickness. 

4.5.  MULTIMODULE  ASSEMBLY 


The  two-  and  three-module  batteries  are  a.ssembled  by  using  the  same  technique 
as  for  a  single  module,  i.e.,  endplate,  insulator,  and  collector  plates  are  positioned 
and  the  cell  stack  is  built  up  in  sequence:  anode,  separator,  cathode  and  bipolar 
plate  -  until  the  required  number  of  cells  are  in  place.  The  first  module  is  completed 
by  placing  a  collector  plate  on  the  last  cathode.  The  second  module  is  built  with 
the  order  of  active  cell  components  reversed,  i.e.,  following  the  sequence:  cathode, 
separator,  anode  and  bipolar  plate.  The  inverted  stacking  order  places  the  high- 
voltage  end  of  one  module  adjacent  to  the  high  voltage  end  of  the  next  module 
and  the  ground  end  of  one  module  to  the  ground  end  of  the  next  module.  This 
arrangement  avoids  a  voltage  gradient  of  the  magnitude  of  several  hundreds  of  volts 
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across  the  collector  plate  when  the  battery  is  activated.  Of  course,  such  high  voltages 
would  cause  extremely  high  intermodule  current  flow  which,  in  all  probability,  would 
lead  to  catastrophic  failure. 


Quality  control  is  normally  achieved  by  measuring  the  impedance  and  phase  an¬ 
gle  of  the  cell  as  a  function  of  frequency.  Well-made  cells  show  capacitances  that 
fall  within  a  standard  deviation  of  5  percent  of  the  mean.  Deviations  greater  than 
this  are  associated  with  a  problem  within  the  cell.  The  method  does  not  permit  a 
determination  of  the  exact  nature  of  the  problem.  However,  the  two  most  common 
causes  are  nonuniform  cell  spacing  and  the  presence  of  internal  shorts. 


5.  CONCLUDING  REMARKS 


In  organising  this  presentation,  we  attempted  to  illustrate  advantages  arising  from 
a  rational  transition  process  when  developing  an  electrochemical  power  source  de¬ 
signed  to  operate  at  its  limits  by  summarizing  our  approach  and  results  based  on 
such  a  transition  philosophy.  Although  we  selected  the  Li/SOC^  couple  for  this 
purpose,  s'.ch  an  approach  is,  undoubtedly,  applicable  to  other  systems  as  well  and 
is  an  example  of  current  trends  in  the  development  of  electrochemical  cells.  In  fact, 
the  discusbion  of  the  concept  of  transition  from  research  to  design  and  construction  of 
electrochemical  devices  was  the  subject  of  a  recent  Palladium  Medal  Address  (New¬ 
man,  1991  I. 


The  thionyl  chloride  electrolyte  system,  both  neutral  and  acidic,  has  a  complex 
compositi  'H  that  determines  the  elementary  reactions  that  occur  on  discharge.  Elec- 
troreducti  m  rarely  involves  free  thionyl  chloride,  rather,  adducts  with  AlCh  and  Li'^ 
undergo  el  xtroreduction  preferential!}'.  The  various  species  involved  in  the  electrore¬ 
duction,  1  )th  reactants  and  products  are  now  known.  However,  the  exact  way  in 
which  the  inal  reactants  fal!  apart  to  form  products,  is  not  yet  clear. 


The  Li/SOCl2  system  has  been  modeled,  quite  successfuly  on  several  levels,  in¬ 
cluding  cell,  module,  and  a.ssembly  of  modules.  Modeling  has  contributed  significantly 
to  an  understanding  of  battery  operation  by  providing  the  basis  for  both  the  design 
and  understanding  of  electrode  thickness  and  porosity,  catalytic  effects,  thermal  man¬ 
agement,  and  intercell  currents.  In  addition,  modeling  guided  the  development  of 
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quality  control  (QC)  procedures  and  considerations. 


QC  is  primarily  a  matter  of  careful  selection  of  materials,  fabrication  of  com¬ 
ponents,  and  assembly  of  these  components  to  the  final  cell,  module,  or  battery 
configuration.  Most  unusual  battery  occurrences  can  be  traced  to  failure  of  one  or 
more  QC  procedures  rather  than,  as  believed  early  in  the  development  cycle,  to  un¬ 
known  chemical  p:  cesses.  As  stated,  modeling  provides  a  reasonable  description 
of  the  characteristics  and  behavior  of  modules  and  batteries.  Specific  designs  are 
proprietary  data  or  are  patent-protected  and  thus  cannot  be  presented. 


LifSOCl2  technology  has  come  a  long  way  since  the  initial  recognition  of  its  value 
as  a  battery  chemistry.  We  believe  that  much  of  this  progress  has  resulted  from  a 
focused  effort  to  pursue  a  logical  progression  of  transitions  through  the  R&D  cycle. 
These  efforts,  as  summarized  above,  have  elucidated  the  chemistry  of  the  system,  the 
optimal  design  features  of  cells  and  modules,  the  assembly  of  these  cells  or  modules  to 
batteries,  and  the  resultant  characteristics  of  these  assemblies.  Again,  Li/SOCh  is 
used  here  only  as  an  example.  Other  battery  technology  developments  could  benefit 
from  such  an  approach. 


6.  LIST  OF  SYMBOLS 


Multiple  use  of  the  same  symbol  for  different  quantities  and  vice  versa  was  em¬ 
ployed  to  conform  with  their  usage  in  quoted  references. 

a  area,  m“^ 

c  concentration,  mole  dm~^ 

Cp  specific  heat,  J  g“^  K"* 

e~  electron,  charge 

D  diffusion  coefficient,  m  s~^ 

E  energy  (activation  energy  in  eq.  (21)) 

F  Faraday  constant,  96487  C/equiv 

AG  Gibbs  free  energy,  J 
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enthalpy,  J 

h  heat  transfer  coefficient,  W  cm~^  K“^ 

I  current,  A 

i  current,  A 

i  running  index 

J  cell  current,  A 

j  current  density,  A  m"^ 

jo  exchange  current  density,  A  m— 2 

K3  solubility  product,  mol^ 

k  rate  constant  (dimension  consistent  with  reaction  order),  eqs.  12,  13, 

14 

/  length  (coordinate),  m 

m  ,  n  running  indices 

m  morphology  parameter,  eq.  13 

p{l,t)  morphology  parameter,  eqs.  11,15 

p  pressure,  atm 

q,  Q  quantity  of  heat,  J 

R  gas  constant,  8.3143  J  mol“*  deg~^ 

R  resistance.  Ohm,  eqs.  25,  26 

r  rate,  s~' 

r  radius,  m 

A5  entropy,  J  deg“* 

T  temperature,  K 

s  specific  surface,  m“',  eqs.  11,  15 

t  time,  s 

U  potential,  V 

V  volume,  m^ 

V  potential,  V,  eqs.  25,  30 

V  volume  fraction,  eq.  11 

V  scan  rate,  mV  s— 1 


68 


W  work,  J 

X  coordinate 

valence 

a  transfer  coefficient,  eq.  10 

0  parameter,  eq.  20 

c  efficiency,  eqs.  11,  15;  sec.  1.3 

C  independent  variable  defined  in  text 

A  independent  variable  defined  in  text,  sec.  2.6 

A  volume  conductivity,  Ohm“^m~',  eq.  18 

A  thermal  conductance,  W  m“^K“^,  eq.  23 

A  parameter,  sec.  2.6;  design  parameter,  sec.  3.5 

Tj  overpotential,  V,  eq.  10 

6  surface  coverage 

i  degree  of  advancement,  eq.  27 

i  parameter,  sec.  2.6;  design  parameter,  sec.  3.5 

cr  source  term  (heat),  J,  eq.  17 

«  thermal  conductivity,  J  cm"*  s"*  deg"*,  eq.  19 

$  potential,  V,  eq.  26 

n  dimensionless  power,  sec.  3.5 

p  dimensionless  resistance,  sec.  3.5 
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APPENDIX  A  -  EXPERIMENTAL 


METHODS  OF  ANALYSIS  OF  THE  BULK  ELECTROLYTE 

Introduction 

Vibrational  spectroscopy  (Raman  and  infrared)  has  be  n  shown  to  be  a  sensitive 
technique  for  probing  solvation  and  association  phenomena  (Popovych  &  Tomkins, 
1981).  Information  on  the  intermolecular  forces,  the  structure  of  solvents,  and  the 
nature  of  solvation  can  be  derived  using  vibrational  spectroscopy.  Solute-solvent  in¬ 
teractions  not  only  affect  the  vibrational  spectra  of  the  solvents,  but  also  lead  to  the 
appearance  of  new  vibrational  bands  due  to  solute-solvent  complexation.  Information 
on  bond  properties  as  well  as  inter  and  intramolecular  interactions  can  be  assessed 
from  the  vibrational  frequencies  and  their  shifts,  band  intensities,  depolarization  ra¬ 
tios,  and  excitation  profiles. 

Instrumentation 


Infrared  spectra  were  obtained  on  a  Nicolet  5DXB  FT-IR  spectrometer  with  a 
resolution  of  4  cm“*.  Barnes  demountable  cells  with  AgCl  windows  and  a  15-;im 
teflon  spacer  were  used. 

The  Raman  spectra  were  recorded  on  a  system  composed  of  a  Lexel  model  85  Ar- 
ion  laser  emitting  at  488  nm,  which  was  chopped  at  20  Hz,  a  sample  chamber  with 
coupling  and  collecting  optics,  a  Spex  model  1400-11  scanning  double  monochroma¬ 
tor  with  photon-counting  preamplifier  and  line  driver,  a  photon-counting  integrator 
(Miller,  1968),  and  a  digital  synchronous  detector.  The  samples  were  contained  in 
a  Helma  model  162  F  stoppered  quartz  cell.  The  cell  chamber  was  flushed  with  dry 
nitrogen  to  prevent  condensation  on  the  cell  windows  when  operating  below  room 
temperature.  Temperatures  were  controlled  to  within  ±0.01°C.  Most  of  the  spectra 
were  recordc;'  at  a  laser  power  of  50  mW. 

Resolution  of  Spectral  Bands 


Raman  and  infrared  vibrational  bands  have  both  a  Lorentzian  and  a  Gaussian 
contribution.  While  the  width  of  the  Lorentzian  line  is  an  intrinsic  property  of  the 
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molecule,  depending  on  the  ratio  of  the  decay  time  to  the  lifetime  of  the  initial  state, 
the  Gaussian  distribution  is  related  to  the  interaction  with  the  environment,  par¬ 
ticularly  with  the  nearest  neighbors  of  the  excited  molecule  (Gabriel,  Mosier-Boss, 
&  Szpak,  1987).  The  observed  excitation  profiles  arise  from  the  convolution  of  the 
Lorentzian  and  Gaussian  contributions  and  assume  the  form  of  the  Voigt  profile 

=  (Al) 

where  a  =  u;  =  2(ir  —  ^Vg  and  no,  Ant  and  Anc 

are,  respectively,  the  center  frequency  and  the  full  widths  at  half  maximum  of  the 
Lorentzian  and  Gaussian  distributioi  >  that  characterize  the  line  (diBartolo,  1968). 


In  a  dynamic  system  composed  of  N  number  of  species,  the  experimentally  ob¬ 
served  spectral  band,  S{u)^  is  a  superposition  of  N  lines  that  are  approximated  by 
the  following  expression 


I{i>)  =  Co{v  -q)  +  I/O,  Ai/Li,  Auoi)  (A2) 

where  the  function  /(...)  is  defined  in  eq.  (Al).  Because  /(...)  is  normalized,  Ci  is  the 
area  of  the  fundamental  Lorentzian  line  of  the  i-th  species.  The  linear  term  Co(i/  —  a) 
accounts  for  a  sloping  baseline. 


The  least  square  fit  of  eq.  (A2)  to  the  experimental  data  is  constructed  by  mini¬ 
mizing  the  test  function,  F,  given  by 

F  =  E",lS(-'r)  -  /(■',)1V(V  -  P)  (A3) 

where  M  is  the  number  of  data  points  in  the  measured  spectrum,  and  P  is 

the  number  of  variable  parameters  characterizing  the  spectral  intensity  calculated  by 
eq.  (A2).  The  minimization  of  F,  eq.  (A3),  was  carried  out  using  the  Nelder-Mead 
algorithm  (Nash,  1979).  The  number  of  lines  N  is  selected  from  consideration  of 
the  vibrational  mode  under  investigation,  the  chemical  nature  of  the  components,  as 
well  as  from  the  minimization  procedure.  Constraining  the  line  widths  to  be  greater 
than  the  data  point  separation  and  including  the  factor  {M  —  P)  in  eq.  (A3),  places 
an  upper  limit  on  the  value  of  N  that  will  minimize  F.  The  integral  appearing  in 
eq.  (Al)  was  evaluated  by  the  method  of  Romberg  quadrature  (Barakat,  1980)  where 


the  4-00  and  — oo  limits  were  replaced  by  4-3  and  —3.  Computational  round  off  errors, 
truncation  errors  in  evaluation  of  eq.  (Al),  and  scatter  in  the  intensity  data  broaden 
the  minimum  in  occasionally  making  it  difficult  to  select  the  global  minimum  from 
several  local  minima.  Except  for  the  truncation  error,  this  difficulty  is  not  specific  to 
the  Voigt  profile,  nor  to  the  Nelder-Mead  algorithm,  but  is  characteristic  of  nonlinear 
optimization. 

Molecular  Orbital  Calculations 

The  MO  calculations  were  performed  using  either  AMPAC  or  MOPAC.  AMPAC 
was  developed  at  the  University  of  Texas,  Austin,  Texas.  MOPAC  is  an  upgraded  ver¬ 
sion  of  AMPAC  developed  by  Dr.  Jimmy  Stewart  at  the  U.S.  Air  Force  Academy,  Col¬ 
orado  Springs,  Colorado.  Both  AMPAC  and  MOPAC  are  general  purpose,  semiem- 
pirical  molecular  orbital  packages.  These  calculations  yield  information  on  the  size, 
shape,  and  charge  distribution  within  molecules  and  complexes. 

METHODS  EMPLOYED  TO  EXAMINE  SOCh 
ELECTROREDUCTION 

In  Situ  Infrared  Spectroscopy 

Electrochemical  Cell.  The  spectroelectrochemical  cell  is  shown  in  figure  Al.  The 
body  of  the  cell  was  made  of  a  machinable  ceramic  material  chemically  resistant  to 
SOCI2  (Macor,  Dow  Chemical).  A  demountable  AgCl  window  was  sealed  against 
the  cell  body  by  an  0-ring  (Kalrez.  DuPont  Chemical).  The  working  electrode  w-as  a 
polished  metal  foil  sealed  to  a  glass  tube  using  a  chemically  resistant  epoxy  (Epoxy 
Patch  C,  Hysol  Corp.).  Ohmic  contact  to  the  metal  foil  was  achieved  using  a  con¬ 
ductive  epoxy.  The  working  electrode  was  held  in  place  in  the  cell  body  by  an  O-riiig 
(Kalrez,  DuPont  Chemical)  and  a  teflon  bushing.  The  counter  and  reference  elec¬ 
trodes  were  a  Pt  wire  ring  and  a  Ag/AgCl  wire  in  a  teflon  sleeve,  respectively.  These 
two  electrodes  were  sealed  to  the  cell  body  (Epoxy  Patch  C.  Hysol  Corp.).  A  variable 
distance  between  the  window  and  working  electrode  surface  was  provided:  however, 
each  set  of  measurements  was  carried  out  at  a  fixed  distance.  The  pathlengths  were 
calibrated  with  neat  CC'U. 


Ag/AgCI 

reference 

electrode 


Working 

electrode 


I  / 

window 


R 

counter 

electrode 


Figure  Al.  Design  of  spectroelectrochemical  cell  used  in  in  situ  IR  reflectance  spectroscopy. 


Instrumentation.  All  in  situ  infrared  spectra  were  obtained  using  the  Nicolet  5DXB 
FT  IR  spectrometer  described  earlier.  The  spectroelectrochemical  cell  was  placed  at 
a  45°  angle  inside  the  spectrometer  using  a  specular  reflectance  apparatus  (Spectra- 
Tech,  Model  500).  Overpotentials  were  applied  using  a  PAR  model  173  potentiostat. 

Relaxation  Techniques 

Electrochemical  Cell.  Figure  A2  shows  the  two-chamber  electrochemical  cell  used 
in  the  relaxation  experiments.  The  body  of  the  cell  is  constructed  of  glass  and  has 
three  glass-threaded  joints  (Ace  Glass).  The  two  chambers  are  separated  by  a  coarse 
glass  frit.  One  chamber  contains  the  working  and  reference  electrodes  while  the 
other  has  the  counter  electrode.  The  working  electrode  was  a  metal  foil  permanently 
mounted  on  a  glass  tube  using  chemically  resistant  epoxy  (Epoxy  Patch  C,  Hysol 
corp.).  Ohmic  contact  was  provided  using  conductive  epoxy.  The  reference  electrode 


was  a  Ag /AgCl  wire  in  a  teflon  sleeve.  Prior  to  the  experiment,  the  reference  electrode 
was  immersed  in  a  fritted  eyedropper  containing  the  electrolyte.  This  prevented  the 
corrosion  of  the  reference  electrode  during  the  course  of  the  experiment.  A  coiled  Pt 
wire  served  as  the  counter  electrode.  The  electrodes  were  positioned  in  the  cell  body 
using  0-rings  (Kalrez,  DuPont  Chemical)  and  teflon  bushings. 
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Figure  A2.  Design  of  electrochemical  cell  used  in  cyclic  voltammetry  and  galvanostatic 
and  potentiostatic  pulsing  experiments. 


Instrumentation.  The  cyclic  voltammograms  and  charging  and  relaxation  curves 
for  the  galvanostatic  and  potentiostatic  pulsing  experiments  were  obtained  with  the 
aid  of  the  computer-driven  potentiostat,  PAR  model  173  with  a  276  IEEE  computer 
interface. 

MODELING  OF  INTERCELL  CURRENTS 


Figure  A3  shows  the  flowchart  of  the  solution  program  u.sed  to  calculate  the  in¬ 
tercell  currents  within  an  assembly  of  modules.  In  this  approach,  the  geometry  and 


physical  properties  of  the  module(s)  are  first  defined.  The  system  parameters  include 
the  number  of  modules  (maximum  of  5),  the  number  of  cells  per  module  (maximum 
of  90),  porosity  fineness  of  the  cathodes,  cathode  use  factor,  anode  reversed  current 
factor  and  load  resistance.  The  cell  and  feed  tube  porosities  are  assigned  as  described 
in  section  3.6.,  using  a  random  number  generator.  The  porosities  were  chosen  from 
the  range  0.3  to  0.9  corresponding  to  limitations  imposed  by  heat  evolution  (Parnell 
&  Szpak,  1985)  and  practical  considerations  of  cell  construction.  In  the  case  of  de¬ 
ficient  modules,  the  number,  location,  and  deficiency  factor  of  the  affected  cells  are 
defined.  In  these  calculations,  the  effective  electrode  spacing  is  assumed  to  be  the 
same  for  each  cell  within  the  module. 


Once  the  system  parameters  have  been  established,  the  starting  cell  and  feed 
tube  conductivities  are  calculated  and  the  time-step  loop  calculations  are  started. 
Within  the  time  loop,  the  port  input  resistances  are  calculated  (Szpak,  Gabriel,  & 
Driscoll,  1984)  and  the  banded  matrix  for  each  module’s  loop  currents  are  set-up.  The 
module  coupling  equations  for  the  coupling  currents  are  computed  from  tridiagonal 
sets  of  linear  equations  -  one  set  for  each  module.  The  reversed  current  for  each 
cell  is  then  computed  as  is  the  charge  transferred  by  the  reversed  current  (Szpak  et 
al.,  1984).  The  loop  currents  are  used  to  calculate  the  charge  transferred  to  each 
cathode.  Using  the  charges  transferred  during  the  time  increment,  the  new  cell  and 
feed  tube  conductivities  are  calculated  (Szpak  et  al.,  1990).  The  time  is  advanced  by 
the  incremental  amount  and  the  time-loop  calculations  are  done  from  the  new  time. 


Figure  A3.  A  flowchart  calculating  the  characteristics  of  the  multiinodule  battery 
assembly,  (Gabriel  et  al.,  1990). 
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these  assemblies.  The  Li/SOClg  is  used  only  as  an  example.  Other  battery  technology  developments  could  benefit  from  such  an 
approach. 
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